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§ 1. Термодинамические характеристики 
химической реакции 

Любая химическая реакция сопровождается выделением 
или поглощением тепла. Количество тепла, которое выде¬ 
ляется или поглощается в ходе химической реакции при за¬ 
данной температуре исходных веществ и продуктов реакции, 
называется тепловым эффектом. 

Тепловой эффект химической реакции, протекающей при 
постоянном объеме, равен изменению внутренней энергии 
системы 

(3 ѵ — АУ- 

Если химическая реакция протекает при постоянном дав¬ 
лении, ее тепловой эффект равен изменению энтальпии 
системы 

<э р = д Н. 

Возможность самопроизвольного протекания реакции 
в заданных условиях определяется знаком изменения соот¬ 
ветствующего термодинамического потенциала: для реакций, 
протекающих при постоянных объеме и температуре, — изо- 
хорного потенциала АР, а при постоянстве давления и темпе¬ 
ратуры — изобарного потенциала АС. 

При V и 7’ = сопз1: Д/ г <0 — самопроизвольный процесс; 
АР = 0 — состояние равновесия; АР > 0 — несамопроизволь¬ 
ный процесс. 

При р и Т = сопзі: ДО <0 — самопроизвольный процесс; 
АО = 0 — состояние равновесия; ДО>0 несамопроизвольный 
процесс. 

Изменение изохорного или изобарного потенциалов равно 
разности между соответствующим тепловым эффектом и 
изменением «связанной» энергии системы: 

АР = Ш— ГД5; 

ДО= АН — ТА5, 

где А5 — изменение энтропии при химической реакции. 
По знаку изменения энтропии при реакции судят о ее на- 
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правлении, но рассматриваемая система в этом случае долж¬ 
на быть изолированной, т. е. должны быть постоянными 
ее объем и внутренняя энергия. Энтропия изолированной 
системы при протекании в ней самопроизвольных процессов 
увеличивается и становится максимальной при достижении 
равновесного состояния. 

Запишем уравнение химической реакции в общем виде 

ѵ Иі ѵ 2^2 ••• ~Ь == Ѵ-Д/ ~Ь ^2 А* + ••• + > (О 

где Аі — химический символ участников реакции; ѵі — соот¬ 
ветствующий стехиометрический коэффициент. Тогда изме¬ 
нение какой-либо из термодинамических характеристик при 
реакции может быть рассчитано с помощью уравнения об¬ 
щего типа 

і і 

где ДУ—изменение термодинамической функции при реак¬ 
ции; У—мольное значение термодинамической функции для 
вещества, участвующего в реакции. 

Более сокращенно последнее уравнение можно записать 
так: 

Д У - 2 ( ѵ Юкон — 2 ( ѵ Юнач, (2) 

где индексы: «кон» и «нач» определяют принадлежность 
участника реакции к продуктам или к исходным веществам. 

По справочным данным о термодинамических свойствах 
индивидуальных веществ рассчитывают основные термодина¬ 
мические характеристики (АН, АЗ и АС) реакции в стандарт¬ 
ных условиях: температура 25°С (298, 15К); давление 
1 атм. 

Для пересчета термодинамической функции с одной тем¬ 
пературы на другую необходимо знать теплоемкости веществ, 
участвующих в реакции. Температурная зависимость тепло¬ 
вого эффекта и изменения энтропии реакции определяются 
алгебраической суммой теплоемкостей, рассчитанной по урав¬ 
нению (2), 

а АН (1&8 Д с„ 

-— Д с • _ —__ 

ат р' ат т • 


Температурная зависимость изменения изобарного 
циала реакции имеет более сложный характер: 


яГДО 

ат 


= -Д5; 


потен- 


причем АЗ представляет собой функцию температуры. 
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Каждая химическая реакция через некоторое время после 
ее начала приходит в равновесное состояние. В этом со¬ 
стоянии состав системы (концентрации или парциальные дав¬ 
ления исходных веществ и продуктов) со временем не изме¬ 
няется. Равновесное состояние химической реакции характе¬ 
ризуется константой равновесия. Наибольшее практическое 
значение имеют константы равновесия, выраженные через 
парциальные давления К ѵ ил-и через концентрации К с - 

Для реакции (1) в соответствии с уравнением изотермы 
реакции 


ДО т 


ДЛп 


Ра/ р а 2 ' 


или 


ДО т 


РаРа-" 

— ЯТ\п К р + ПТ \п 



г ѵ ' ѵ ' 

рі’РЬ-" 

+ РТ1 п 

л \ 2 

ѵ . ѵ_ 

_ равн 

1 

"Ч* 


Ра^РЛ’--- 

РаРа-- 


Если все вещества, участвующие в реакции (1)—исходные 
и конечные, находятся в стандартных состояниях, то в этом 
случае 


АО т = АС т °^-КТ\пК р , (3) 


где АСт° — стандартное изменение изобарного потенциала 
или, сокращенно, стандартный изобарный потенциал реакции. 
Из уравнения (3) 


1п К р = 


ДО т о 
КТ 1 


Соответственно 


1п К с — — 


Д/^о 

КТ * 


В идеальных условиях константы Кр и К с связаны между 
собой простыми соотношениями: 

К р = К с (КТ ) 4ѵ и К с — К р (КТ )~ Лѵ , 

где Аѵ — изменение числа молей газообразных веществ 
в соответствии со стехиометрическим уравнением реакции. 
Для реакции (1), если все вещества являются газообраз¬ 
ными, 

д/ аз = 2Ѵ-2ѵ 

В справочниках, как правило, приводятся данные, позво¬ 
ляющие рассчитывать термодинамические характеристики 
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реакции в изобарно-изотермических условиях ( р, Г=сопзі). 
В некоторых случаях может возникнуть необходимость рас¬ 
чета термодинамики реакции в изохорно-изотермических 
условиях (например, синтез полупроводниковых соединений 
часто проводится в запаянных ампулах). В этом случае тер¬ 
модинамические функции АИ и АР, характеризующие проте¬ 
кание реакции в условиях постоянства температуры и объе¬ 
ма, могут быть вычислены из своих термодинамических ана¬ 
логов АН и АО по формулам: 

ДІ7 = АН — 7"Дѵ газ ; 

ДЯ = ДО — /?ГДѵ газ , 



где Дѵ газ — алгебраическая сумма стехиометрических коэф¬ 
фициентов газообразных участников реакции. 

Значение константы равновесия химической реакции поз¬ 
воляет рассчитать количественную сторону химического пре¬ 
вращения, т. е. количество веществ, вступающих в реакцию 
и образующихся в ходе ее протекания. Обычно степень пол¬ 
ноты реакции характеризуют степенью превращения одного 
из реагирующих веществ. 

Степень превращения — это отношение количества веще¬ 
ства, вступившего в реакцию п', к исходному количеству этого 
вещества щ 

п’ 

Как видно из этого определения, степень превращения яв¬ 
ляется своеобразным коэффициентом полезного действия 
реакции. Если химическая реакция заключается в распаде 
какого-нибудь вещества, степень превращения обычно име¬ 
нуют степенью диссоциации. 

Степень превращения зависит как от константы равнове¬ 
сия, так и от соотношения между исходными количествами 
реагентов. Конкретный вид связи между этими величинами 
определяется стехиометрией химической реакции. 

§ 2. Расчет теплового эффекта реакции 

Исходными данными для расчета теплового эффекта хи¬ 
мической реакции служат теплоты образования веществ, 
участвующих в реакции. В справочниках приводятся тепло¬ 
ты образования химических соединений в стандартных усло¬ 
виях, обозначаемые через АН® 298 и используемые для рас¬ 
чета теплового эффекта химической реакции при 298 К по 
уравнению (2) 

^■^293 = 298)кон (^рАН®, 29в) н ач> 
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где ѵ — стехиометрические коэффициенты в уравнении ре¬ 
акции. 

Пересчет теплового эффекта реакции с температуры 298 К 
на любую другую температуру осуществляется по уравнению 
Кирхгофа, которое в интегральной форме имеет вид 

ДЯ т ° = Д//° 8 + ] Ас р (ІТ, 

298 

где 

= (2 ѵ ^р)кон (2 ѵ ^р)нач> 

Ас р — алгебраическая сумма теплоемкостей для данной 
реакции. 

Если внутри рассматриваемого температурного интервала 
один или несколько участников реакции испытывают фазо¬ 
вые превращения, температурный интервал разбивают на не¬ 
сколько частей 

г т 

провр г 

Д//т° = Д//° 98 + | Д с/аТ±АН прпр + } д с” АТ. 

298 ^превр 

В каждую алгебраическую сумму А с/ и А с р " входят тепло¬ 
емкости тех фаз, которые устойчивы в данном температур¬ 
ном интервале. Теплота превращения АЯщевр берется со 
своим или противоположным знаком в зависимости от при¬ 
надлежности к продукту реакции или к исходному веществу. 

*» 

Методы определения теплоты образования 

Как известно, под теплотой образования соединения пони¬ 
мают тепловой эффект реакции образования этого соедине- 
' ния из свободных элементов в стандартных условиях, т. е. 
при давлении 1 атм и температуре 25° С. Отсюда следует, что 
теплота образования элемента в стандартных условиях равна 
нулю. 

Теплоты образования определяются калориметрически, и 
данные о них содержатся в физико-химических справочни¬ 
ках. Если справочные данные для рассматриваемого соеди¬ 
нения отсутствуют, то теплоту образования можно оценить 
с помощью одного из нижеприведенных методов. 

Метод Беркенгейма. В основе метода лежит положение, 
что в одной и той же группе или в одном и том же ряду 
периодической системы теплота образования соединения 
равна среднему арифметическому из теплот образования ана¬ 
логичных соединений элементов, соседних с данным в ряду 
или периоде таблицы Менделеева. Например, Мд стоит 
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между Ыа и А1 в третьем периоде таблицы Менделеева, 
тогда 


ДЯ° 298 (М§С1 2 ) 


298 (N801) + ЬН/, 298 ( АІСІ3) 

: 2 


(-97,9)+ (-53,7) 
2 


= — 75,8 ккал/моль. 


Экспериментальное значение теплоты образования хлори¬ 
стого магния составляет —75,5 ккал/моль. 

Метод Лотье — Карапетьянца. В рядах подобных веществ 
теплоты образования изменяются по линейному закону. На¬ 
пример, если по оси абсцисс 
Ш-\ / отложить теплоты образо¬ 


зоо- 


С 

•ѵ) 



■ т т 

~ Ь-Н * 1 гзе(5Ч'Гі І ) 


Рис. 1. Определение теплоты образо¬ 
вания $іС1 4 аз методом Лотье — Ка¬ 
рапетьянца 


вания галогенидов олова, 
а по оси ординат —- соот¬ 
ветствующие им теплоты 
образования галогенидов 
кремния, то полученные точ¬ 
ки достаточно хорошо укла¬ 
дываются на прямую линию 
(рис. 1). Отсюда вытекает 
следующая методика опре¬ 
деления теплоты образова¬ 
ния. Предположим, нас ин¬ 
тересует теплота образова¬ 
ния газообразного четырех¬ 
хлористого кремния. Пусть 
мы располагаем данными о 
теилотах образования газо¬ 
образных ЗіР 4 и 5іВг 4 . 
Таким образом, имеем сле¬ 
дующий ряд соединений 
5 іЕ 4 — 5іС1 4 — ЗіВг 4 . Теперь 
необходимо отыскать второй 
ряд, ему подобный. Такими 
рядами могут быть: 


СР 4 — СС1 4 — СВг 4 ; ОеР 4 — ОеС1 4 — ОеВг 4 ; 
ЗпР 4 — ЗпС1 4 — ЗпВг 4 . 


Рассмотрим последний ряд. Строим график: по оси абсцисс 
откладываем теплоты образования галогенидов олова, а но 
оси ординат — соответствующие им теплоты образования 
галогенидов кремния (рис. 1). Через две точки, соответствую¬ 
щие теплотам образования фторидов и бромидов, проводим 
прямую линию. Еще лучше было бы привлечь данные и по 
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йодидам, тогда прямую можно было бы провести по трем 
точкам. Теперь из абсциссы, соответствующей теплоте обра¬ 
зования ЗпС1 4 , восстанавливаем перпендикуляр до пересече¬ 
ния с прямой и по ординате отсчитываем теплоту образова¬ 
ния 5 іС 1 4 , получаем— 150 ккал/моль, что хорошо согласуется 
с экспериментальным значением—156,4 ккал/моль. 

Расчет теплоты образования по энергиям связи. Синтез 
соединения из простых веществ можно отождествить с про¬ 
цессом последовательного перевода элементов в одноатомные 
газы и образования из них соединения. Первый этап связан 
с затратой энергии на разрыв связей, а второй — с выделе¬ 
нием энергии образования новых связей. 

Рассмотрим теплоту образования газообразного ЗіС1 4 , 
т. е. тепловой эффект реакции 

5і тв + 2С1™ 3 = ЗіС1 4 аз ДЯ° 298 (5іС1$ аз ). (4) 

С другой стороны, образование ЗіС1 4 можно было бы про¬ 
вести следующим путем: 

а) атомизации кремния 

5і тв = Зі газ ДЯ° 298 (5Р 3 ); (5) 

б) атомизации хлора 

1/2С1 г « = СГ аз Д#° 298 (СГ аз ); (6) 

в) образования 5іС1 4 из атомов 

5і газ + 4СГ аз = 51С1 Г / 3 ДЯ° = - 4 Е (31 - С1). 

Последний тепловой эффект равен энергии образования че¬ 
тырех связей 5і — С1; ей приписывается знак минус, посколь¬ 
ку энергия при этом выделяется (теряется системой). 

Если уравнение (6) умножить на 4 и сложить последние 
три уравнения, то мы придем к исходному уравнению (4). 
Таким образом, по з'акону Гесса 

ДЯО 298 (5ІС1$ аз ) = Д Щ 298 (5і газ ) + 4ДЯ° 298 (СГ 3 ) - 
— 4Е (8і — С1). 

Энергия связи 5і — С1 равна 97 ккал/моль, теплоты атоми¬ 
зации кремния и хлора соответственно равны: 

ДЯ“ 298 (5і газ ) = 112 ккал/моль; 

Д//° 298 (СГ аз ) = 29 ккал/моль. 

Таким образом, теплота образования хлорида кремния 

Д Щ „ 98 (ВіС1 ваз ) =112 + 4|29 — 4|97 = - 160 ккал/моль. 
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Экспериментальное значение 

ДЯ° 298 (31С1-) = — 156,4 ккал/моль, 

т. е. расхождение составляет 4%. 

В общем случае образования любого соединения 

д/у обр = 2м//° 298 (э газ ) - 

I і 

где Е { — энергия связей; Д Н°^ т (Э т&3 ) — теплоты атомизации 
элементов; п { — число атомов; /п, — число связей. 

Теплоты ат омизации элементов приводятся в справочни¬ 
ках. Иногда они именуютс я теплотами образования атомов. 
а иногда теплотами образования одноатомных газов из эле- 


Определение теплот фазовых переходов 

Если справочные данные по теплотам агрегатных превра¬ 
щений отсутствуют, можно воспользоваться одним из ниже¬ 
приведенных приближенных методов. 

Правило Тритона. Интегрирование уравнения Клапейро¬ 
на — Клаузиуса для процесса испарения в приближении по¬ 
стоянства теплоты испарения приводит к выражению 

1п р^-^ + С, 


где р — давление насыщенного пара; С — постоянная инте¬ 
грирования. 

В точке кипения жидкости 

Т Т кпП , р= 1 атм, А// Исп = Д// кип . 
Следовательно, 


Г' _ АЯкщі 
Ьт * 

* КИП 


откуда 


Эта формула известна в литературе под названием правила 
Трутона. Правило имеет приближенный характер, однако для 
очень многих жидкостей 

ѵ о. кал /„о Дж \ 


моль-град 


моль-град 


Отклонения от правила Трутона обнаруживают жидкости 
с низкой температурой кцпения,, а также ассоциированные 
жидкости. 


Для более точной оценки теплоты испарения по правилу 
Трутона рекомендуем нахождение величины К по химиче¬ 
скому аналогу. Например, оценим теплоты испарения в точке 
кипения ВЛ 3 (Г кип = 483 К), если известно, что для его фи¬ 
зико-химического аналога ВВг 3 ДЯ кип — 7,3 ккал/моль при 
Т'кип = 364 К. Для ВВг 3 


7 кип 


7300 

364 


20,055. 


Тогда для ВЛ 3 

Д// кип = 483• 20,055 — 9,7 ккал/моль. 


Экспериментальное значение теплоты испарения ВЛ 3 в точке 
кипения составляет 10,0 ккал/моль (расхождение 3%). 

Указанная оценка стано¬ 
вится еще более надежной, 
если она опирается на дан- 
ные по теплоте испарения и 

нескольких аналогов. Прак¬ 
тически удобнее построить ю 

зависимость ДЯкип от Т кип 
для нескольких подобных 3 

веществ и из нее найти теп- е 

лоту испарения исследуе¬ 
мого вещества. Для приме- 7 

ра на рис. 2 построена за¬ 
висимость ДЯкип от Т кѵт для 6 

нескольких галогенидов 5і, Ш ХО Ы 

Ое И 5п. Найденные ИЗ нее р ис 2. Зависимость теплоты испаре- 
теплоты испарения ОеСІ 4 ния в точке кипения от температуры 
(7,5 ккал/моль) и ОеВг 4 кипения в ряду подобных веществ 
(9,9 ккал/моль) очень хо¬ 

рошо согласуются с экспериментальными величинами (соот¬ 
ветственно 7,5 и 9,7 ккал/моль). 

Теплоты плавления тоже подчиняются сходному правилу, 
однако величина К. остается постоянной только в пределах 
одного класса подобных веществ и меняется при переходе 
от органических вещестб к неорганическим или от простых 
веществ к сложным. Для неорганических соединений можно 
принять 



ДЯ ПЛ ^ „ кал /„с Дж 
Т пл моль-град ѵ моль • град 


По абсолютной величине теплота плавления всегда зна¬ 
чительно меньше теплоты испарения, а потому ошибка в ее 
определении мало сказывается на точности конечного резуль¬ 
тата определения теплового эффекта химической реакции. 
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Оценка теплоты испарения жидко- 
„ _>е может быть произведена по пра¬ 
вилу Гуггенгейма, согласно которому энтропии испарения 
подобных веществ при одинаковой приведенной температуре 
имеют одинаковое численное значение (приведенной темпе¬ 
ратурой называют отношение заданной температуры к кри¬ 
тической температуре вещества). 

Из правила Гуггенгейма следует, что зависимость энтро¬ 
пии испарения от приведенной температуры для подобных 

веществ должна описы- 
д5исп . ваться одной и той же 

кривой. Отсюда для ве¬ 
щества, теплоту испаре¬ 
ния которого нужно опре¬ 
делить, подбирается не¬ 
сколько подобных ве¬ 
ществ, теплота испарения 
которых заранее извест¬ 
на. Обычно значение теп¬ 
лоты испарения известно 
для температуры кипе¬ 
ния. Тогда по закону 
Кирхгофа 

/1\Н иС11 _ п г ж 

ЛТ ~~ с р С Р 

рассчитываются, темпера¬ 
турная зависимость теп¬ 
лоты, а затем и энтропии 
Рис. 3. Зависимость изменения энтропии испарения, поскольку 
при испарении от приведенной темпера- 

туры в ряду подобных веществ А5 ИСП = - сн -. • 

Г исп 

Далее, зная критическую температуру, вычисляют приведен¬ 
ные температуры и строят график 

Д5 И с П =/(7’Д. 

Для вещества, теплоту испарения которого необходимо 
найти, рассчитывают приведенную температуру, исходя 
из того значения температуры, при котором следует оце¬ 
нить Д# исп . По графику Д 8псп = !(Т Г ) для вычисленного 
значения Т г определяют А5 ИСП , а затем пересчитывают ее 
на АЯ ИСП . 

Например, определим теплоту испарения в точке кипе¬ 
ния для ВіСІз. В качестве подобных веществ выбираем РСІз, 
АзС1 3 и 5ЬС1з. Для этих веществ строим зависимость Д5 ИСП ?= 
= І(Т Г ) (рис. 3). Используя Г КШ і = 714К и 7’, Ф =1178К хло- 
12 



Правило Гуггенгейма. 
сти при любой температуг 


рида висмута, подсчитываем приведенную температуру этого 
вещества - , • • ■ 



По графику на рис. 2 для этого значения приведенной темпе¬ 
ратуры находим Д$ исп = 22,5 кал/моль • град. Отсюда тепло¬ 
та испарения в точке кипения ВіС1 3 

Д// исп = 714-22,5 = 16,08 ккал/моль. 

Экспериментальное значение равно 17,35 ккал/моль (рас¬ 
хождение 7,3%). 


§ 3. Расчет изменения энтропии реакции 

Изменение энтропии реакции при стандартных условиях 
подсчитывается как алгебраическая сумма абсолютных 
энтропий участников реакции 



Пересчет изменения энтропии с температуры 298 К на любую 
другую осуществляется по формуле 

т 

М> т =№ ш +\^йТ. 

298 


При наличии фазовых переходов интервал интегрирования 
разбивается на отдельные участки 



Если справочные данные по энтропии веществ отсут¬ 
ствуют, то величину энтропии можно оценить по одному из 
следующих приближений. 

Для твердых тел: 

— по Беннеру энтропия твердых соединений линейно за¬ 
висит аГлогарйфма молекулярного веса 

5 = Л1 ё М + Д, 


где Л и В — постоянные для однотипных веществ. Поэтому, 
имея значения энтропии двух веществ какой-либо группы 
соединений (например, окислов металлов), можно найти не¬ 
известную энтропию других соединений этой группы; 
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— по Герцу для элементов, галоидных солей и сульфидов 
справедливо уравнение 

з г _ 

8 у с р = сопзі — 4,9 п, 

где с р — теплоемкость при 298 К, а п — число атомов в мо¬ 
лекуле соединения. 

Если для данного класса химических соединений отсут¬ 
ствует эмпирическая формула для оценки энтропии, то ее 
приближенно можно определить простым суммированием 
энтропий составляющих элементов. Например, для соедине¬ 
ния Мо 3 Зі: 

5 (Мо) = 6,83 кал/моль • град; 8 (5і) — 4,50 кал/моль • град; 

8 (Мо 3 8і) — 3 • 6,83 + 4,50 = 25,0 кал/моль • град. 

Экспериментальное значение энтропии Мо 3 5і равно 
25,4 кал/моль • град. 

Для газов: 

— для двухатомных газов значение энтропии может быть 
найдено по эмпирическому уравнению 

5 298 = 220,8 + 0,180 М - 1005 М -1 , 
где М — молекулярный вес; 

— для многоатомных газов энтропии оцениваются по 
уравнению 

5 298 = 163 + 1,42 М - 26,0 ■ 10 -4 М 2 . 

Для приближенного определения энтропии жидких ве¬ 
ществ можно воспользоваться энтропией этих веществ в твер¬ 
дом или газообразном состоянии. Тогда 

5 ж = 5 тв + Д5 пл ; 

5 ж сгаз до 

О а *^ИСП) 

где А5 ПЛ и Д5 ИСП — изменение энтропии при плавлении и 
испарении. 

Для веществ в любом фазовом состоянии энтропия может 
быть определена по методу Лотье — Карапетьянца. В рядах 
подобных веществ значения энтропии при стандартных усло¬ 
виях связаны простыми линейными соотношениями: 

№) і = л + 5 ( 5 2° 9 8)п; (7) 

(-^гэв)I — й ^ ( С р , 298)і> 

где индексы I и II относятся к рядам однотипных соедине¬ 
ний; А, В, а и Ь — эмпирические постоянные для данных 
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рядов. Например, пусть необходимо определить 5^ (ЗіСі^* 3 ). 
Выберем ряд однотипных соединений 

ССЦ - ЗіС1 4 - ОеСЦ. 

Для того чтобы найти $® 98 ( ЗіСІ^ 3 ), необходимо знать с® 298 
всех соединений, а также стандартные энтропии СС1 4 и 
ОеСЦ. Тогда можно взять данные для двух веществ, соста¬ 
вить два уравнения типа (7) и найти постоянные а и Ь. 
Более просто и наглядно определение 3^ 98 (ЗіС1 4 ) можно 



Рис. 4. Оценка энтропии 5іСЦ Рис. 5. Оценка энтропии 8іСЦ ме- 
методом Лотье — Карапетьянца тодом Лотье — Карапетьянца 


провести графическим путем. Для этого строим график 
8 = [(Ср) в ряду соединений ССЦ— ЗіС1 4 — ОеСЦ, отклады¬ 
вая по оси абсцисс с® 298 , а по оси ординат соответствующие 
им 5^ 98 (рис. 4). Полученные точки соединяем прямой ли¬ 
нией. Теперь восстанавливаем перпендикуляр из абсциссы, 
соответствующей с® т ( ЗіСІ™ 3 ), до пересечения с прямой 
линией и на оси ординат находим 5^ 98 (ЗіС1 4 аз ). Эту же 
величину можно определить и из рассмотрения изменения 
энтропии в рядах ССЦ— ЗіС1 4 —ОеСЦ и СВг 4 — ЗіВг 4 —ОеВг 4 
(рис. 5). Графическая оценка 5^ 8 (ЗіСІ™^ кал/моль-град) 
по рис. 4 дает 78,7, а по рис. 5 — 78,8, что хорошо согласуется 
с экспериментальным значением 79,1 кал/моль • град. 
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§ 4. Расчет константы равновесия химической реакции 

Главной целью термодинамических расчетов является 
определение изменения изобарного потенциала химической 
реакции. По имеющимся в справочной литературе данным 
можно рассчитать стандартное изменение изобарного потен¬ 
циала АО 0 , относящееся к состоянию реакции, в котором 
парциальное давление каждого газообразного вещества рав¬ 
но 1 атм. Стандартное изменение изобарного потенциала по¬ 
зволяет определить главную характеристику химической реак¬ 
ции-константу химического равновесия по уравнению (3). 
Зная численное значение константы равновесия, можно рас¬ 
считать выход любого продукта реакции и оценить оптималь¬ 
ные условия для проведения реакции. 

Кроме того, зная АС 0 , можно произвести оценку термо¬ 
динамической вероятности протекания реакции. Однако 
к такой оценке следует подходить с осторожностью. Абсо¬ 
лютно достоверной характеристикой возможности протекания 
химической реакции в данных условиях служит истинное 
изменение изобарного потенциала АС, связанное со стандарт¬ 
ным изменением АС 0 уравнением изотермы химической реак¬ 
ции (стр. 5). В общем случае знак АС может не совпадать 
со знаком АС 0 , а потому последний является лишь прибли¬ 
женным критерием возможности протекания реакции. Однако 
правильные выводы о термодинамической возможности реак¬ 
ции можно сделать и по знаку АС 0 , если давления в рассмат¬ 
риваемой системе не слишком сильно отличаются от 1 атм, 
или абсолютное значение АС 0 достаточно велико (10 и более 
ккал или 40 и более кДж) . Существуют различные методы 
расчета стандартного изменения изобарного потенциала реак¬ 
ции, из которых наиболее употребительными являются сле¬ 
дующие. 

Метод комбинирования химических равновесий 

Сущность метода заключается в подборе двух или не¬ 
скольких уравнений химических реакций с известными термо¬ 
динамическими характеристиками, сложение которых могло 
бы дать уравнение исследуемой химической реакции. 

Поскольку изобарный потенциал является функцией со¬ 
стояния системы, т. е. его изменение определяется только 
начальным состоянием (состояние исходных реагирующих 
веществ) и конечным состоянием (состояние продуктов реак¬ 
ции) системы и не зависит от пути перехода, численное зна¬ 
чение АО 0 может быть определено как результат соответ¬ 
ствующих алгебраических вычислений. 
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В качестве примера рассмотрим определение изменения 
изобарного потенциала реакции 

5іСі;“ + 2Н г 2 аэ = 5і тв + 4НСГ аз ДО 0 . 

Уравнение этой реакции можно получить из следующих двух 
реакций: 

5Г В + 2С1” 3 = ЗіС1^ аз ДО,; 

Н$ аз + СІ5 аз = 2НСГ аз Д0 2 . 

Однако для этого первую реакцию следует домножить на —!, 
а вторую — на + 2 и после этого сложить. Соответственно 
АО 0 суммарной реакции: 

ДО 0 — - ДО, 0 -4-2ДО 2 0 . 

При вычислении константы равновесия суммарной реак¬ 
ции по константам равновесия составляющих реакций необ¬ 
ходимо помнить, что АС 0 связано с К Р логарифмической за¬ 
висимостью, а потому умножение А 0° соответствует возведе¬ 
нию Кр в степень сомножителя. Тогда 


Энтропийный метод расчета 

После того, как на основании постулата Планка были 
определены абсолютные значения энтропий и составлены 
сводки стандартных значений этой величины для многих ве¬ 
ществ, возможен расчет АС 0 по стандартным теплотам обра¬ 
зования и стандартным энтропиям веществ, участвующих 
в.реакции. В основе расчета лежит известное термодинами¬ 
ческое уравнение , 

ДО 0 =ДЯ°- 7*5°.. 

Исходными данными для расчета служат стандартные теп¬ 
лоты образования Д//^ 298 и стандартные энтропии ве¬ 
ществ 5 2 ° 98 . 

На первом этапе определяют изменение энтальпии реак¬ 
ции при 298 К 

Д ^298 = (ЕѵДТ/у? 298) кон “ 298)нач 

и изменение энтропии реакции при этой же температуре 


Д5 2 ° 98 =№) 


^29в)нач' 


2 
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В этих уравнениях ѵ — стехиометрические коэффициенты; 
Д/^_ 298 и 3 ж —мольные значения соответствующих термо¬ 
динамических функций. Полученные значения Д#^ 98 и Д5» 98 
пересчитывают на требуемую температуру с помощью алгеб¬ 
раической суммы теплоемкостей веществ, участвующих в хи¬ 
мической реакции 

~ (^ ѵ ^)кон (^ Ѵ( -р)нач> 

по формулам: 

д//»=дя» 8 +Г д С р ат ; 

298 

Т 

Д5» = Д5 2 ° 98 +1 ~ ат. 

29* 

Если в исследуемом интервале температур участники 
реакции претерпевают фазовые превращения, то интервал 
интегрирования разбивается на отдельные участки: 

Г прсвр I т 

Д#»=Д //» 8 + | кСр'йТ ± ДЯ преір + _) Д с/ат; 

298 Г превр 

Г прсвр г Т 

Д5» = Д5» + Г -^1 ат+ - А -.{ П -Р с “ |) + ^р’іат. 

•- т - т пвев0 ‘ ^ г 


Теплота превращения берется со своим знаком, если она 
относится к продукту реакции, и с противоположным, если 
она относится к исходному веществу. 

Расчет изобарного потенциала можно провести в трех 
приближениях. 

В первом, самом грубом, приближении пренебрегают тем¬ 
пературной зависимостью термодинамических функций, т. е. 
полагают Дс р ^=0. Тогда 

ДО 0 = Д//° — ГД*? 0 

и// 298 1 и °298 - 

При втором приближении полагают 
Д^ = сопзі = Дс /) 298 . 


Тогда 


- ^^ 98 + Дс р , 298 ( Т - 298 - Лп 2^) 
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Второе приближение энтропийного метода в литературе 
известно под названием приближения Улиха, который придал 
последней формуле вид 

40? - - гая?» - и ,.жТ 0" ш + ™ - 1 ) 

и, обозначив сомножитель 

/ (Т) ІП 298 і 1 > 

рассчитал его численное значение для разных температур. 

Таблица 1 


Функция Улиха Т/(Т) — И——^— — 1^ 


т, к 

Т/(Т) 

Т, К 

ЩТ) 

300 

0,00 

1700 

1557,48 

400 

15,69 

1800 

1734,45 

500 

56,65 

1900 

1916,97 

600 

117,74 

2000 

2104,76 

700 

195,58 

2100 

2297,55 

800 

287,76 

2200 

2495,10 

900 

392,46 

2300 

2697,20 

1000 

508,30 

2400 

2903,65 

1100 

634,16 

2500 

3114,27 

1200 

769,12 

2600 

3328,89 

1300 

912,42 

2700 

3547,35 

1400 

1063,43 

2800 

3769,52 

1500 

1221,58 

2900 

3995,27 

1600 

1386,40 

3000 

4224,46 


Значения }{Т) и произведения Т ((Т) сведены Улихом 
в табл. 1. Таким образом, расчет изобарного потенциала по 
методу Улиха ведется по формуле 

ДО» = ДЯ°, 8 - ГД5 298 - д с р , ш т/ (Г). 

При третьем,, наиболее точном, приближении учитывают 
температурную зависимость термодинамических функций. 
В качестве примера рассмотрим расчет термодинамики 
реакции 

ЗіС15 аз + 2Щ аз = Зі тв + 4НСГ аз (8) 

энтропийным способом для Г= 1000К. 

Из справочников находим основные свойства веществ, 
участвующих в реакции (8) (табл. 2). 

2 * 
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Таблица 2 


Термодинамические характеристики реагентов 


Вещество 

298 * 

ккал 

с0 

° 298» 

кал 

с р = « + РТ - 

к ал /мо ль -г 

і 

, т— 2 

7 Т 
■рад 

С р> 298» 
кал 


моль 

моль -град 

а 

р-103 

і' ■ 10~ 5 

моль• град 

81С17 3 

-156,4 

79,10 

24,25 

1,64 

2,75 

21,65 

ілгаз 

Н 2 

о 

31,21 

6,62 

0,81 

— 

7,06 

ЗЦВ 

0 

4,50 

5,66 

0,73 

1,04 

4,71 

НС1™ 3 

-22,1 

44,62 

6,70 

0,84 

— 

6,95 


Тепловой эффект реакции при 298 К 

ДЯ° 8 = 4ДЯ° 298 (НСГ 3 )- Д//« 298 (8іС1- 3 ) =4(—22,1) - 

— (— 156,4) = 68,0 ккал. 

Изменение энтропии реакции при 298 К 
Д^298 = 45 2 ° 98 (НСР 3 ) + 5 2 ° 98 (5і тв ) - 5 2 ° 98 (8ІС1-) - 25« 98 (Н-) = 

= 4-44,62 + 4,50 -79,10 -2-31,21 = 182,98-141,52 = 

= 41,46 кал/град. 

Первое приближение. Полагаем Дс р = 0, тогда 
\И° — ДЯ° • Д < ?° = Д < ?° • дГ?° = ДЯ 0 тдсо • 

а/і 298> Т а °298’ Т а/і 298 1 ао 298» 

АО шю = 68 000 - 1000 • 41,46 = 26 540 кал. 

Второе приближение. Для расчета во втором приближе¬ 
нии необходимо подсчитать алгебраическую сумму теплоем¬ 
костей для температуры 298 К 

Д*р,яи = 4,71 +4-6,95-21,65-2-7,06 = 32,51 - 35,77 = 

= — 3,26 кал/моль-град. 

Тогда 

Д Я» 000 = ДЯ° 98 + Д с р> 298 (1000- 298) = 68 000 - 3,26 • 702 = 

= 65 710 кал; Д5 1000 =Д5 2 ° 98 + Дс р>298 1п-^- = 

= 41,46 - 3,261п -Ш- = 41,46 - 3,95 = 37,51 кал/град. 
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Таким образом, во втором приближении 

ЛО і°ооо = 65710 - 1000-37,51;= 28 200 кал. 

Третье приближение. Рассчитываем алгебраическую сумму 
теплоемкостей как функцию температуры, используя для 
каждого вещества соответствующий температурный ряд: 

А+ = — 5,03 -ф-0,83• 10 -3 Г + 1,71 • 10 5 7'~ 2 кал/моль-град. 
Тогда 

1000 

дя°оо = дя°8+ | д^г= 

1000 298 

= 68000+1 (—5,03 + 0,83-10~ 3 7' + 1,7М0 б Г -2 )=65250 кал; 

298 

1000 

Д^ооо = Д+9 8 + | %-ат= Д5 2 ° 98 + 


+ 1 


— 5,03 + 0 , 83-10 _3 7 + 1,71 • 10 5 7'~ 2 


<ІТ — 36,82 кал/град. 


Следовательно, 

д °і°ооо = 65 250 -1000- 36,82 = 28 430 кал. 


Итак, по мере уточнения значения Дменяется сле¬ 
дующим образом: в первом приближении равно 26 540 кал; 
во втором —28 200 кал, а в третьем — 28 430 кал. 


Метод Темкина и Шварцмана 

Если в рассматриваемом интервале температур вещества, 
участвующие в реакции, не испытывают фазовых переходов, 
расчет энтропийным методом значительно ускоряют, исполь¬ 
зуя специальные вспомогательные таблицы, составленные 
Темкиным и Шварцманом. 

В формуле для изменения изобарного потенциала реакции 

т т 

д 0°т = АЯ° 98 - 7*55* + | а с р ат- т\^-йТ 

298 298 

последние два слагаемых можно объединить в одно, исполь¬ 
зуя формулу интегрирования по частям. Как известно, 

| й («• ц) = | (ийѵ + ѵйи) 

или 


и-ѵ = | ийѵ + ^ѵйи. 


21 


Если положить 


И =І- 


и ѵ = Т, 




,Дс ^ г лт. 


ѵйи =\т 


ьс р ат. 


Таким образом, 


Тогда 




Д6'0 = ДЯ« 8 - 7*5°* + ис р сіт 

298 

- Г Щг Г д С р ат- \ Ас р ат 

298 298 298 


и, следовательно, третий и пятый члены в правой части урав¬ 
нения сокращаются. 

Разделив все члены на Т, будем иметь 


до® міі. 


аа Л СО 

Т "*°298 


т Г 

\ и р лт. 


298 298 


Как видно, последний член правой части уравнения дает по¬ 
правку к результатам приближенного расчета. Введем обозна¬ 
чение 


т т 

г йТ г 


и= | ТЧТ. 


Тогда при п = О 


298 298 


л* і Т . 298 . 

/Кд ІП 2 дд I р 1 , 


а при пф О 


м _ Т п , 298 л+1 298” 

” л(л + 1) "*■ [П+ 1) Т п • 
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Таблица 3 


Величины М п для вычисления термодинамических функций 
по методу Темкина и Шварцмана 

.. , Т 298,15 

М ° - Іо 298>15 + -у— - 1 


г Т п 298,15 я+1 298,15” 

” - Л (п+ 1) + (И + 1 ) т - — Т~ {П Р И п * °) 

^-ю -3 М,.10~ 6 м_ г ІО 5 М г ІО -9 


300 0,0000 0,0000 

400 0,0392 0,0130 

500 0,1133 0,0407 

600 0,1962 0,0759 

700 0,2794 0,1153 

800 0,3597 0,1574 

900 0 4361 0,2012 

1000 ѵ 0,5(588 0,2463 4 

1100 * 0,5765 0,2922 

1200 0,6410 0,3389 

1300 0,7019 0,3860 

1400 0,7595 0,4336 

1500 0,8141 0,4814 

1600 0,8665 0,5296 

1700- 0,9162 0,5780 

1800 Т),9635 0,6265 

, 1900 1,009 0,6752 

2000 1,0525 0,7240 

2100 1,094 0,7730 

2200 1,134 I 0,8220 

2300 1,173 0,8711 

2400 1,210 0,9203 

2500 1,246 0,9696 

2600 1,280 1,0189 

2700 1,314 1,0683 

2800 1,346 1,1177 

2900 1,3775 1,1672 

3000 1,408 1,2166 


0,0000 0,0000 0,00000 

0,0043 0,0364 0,00144 

0,1490 0,0916 0,00553 

0,0303 0,1423 0,01246 

0,0498 0,1853 0,02257 

0,0733 0,2213 0,03630 

0,1004 0,2521 0,05411 

0,1134 0,2783 0,07647 

0,1652 0,2988 0,10388 

0,2029 0,3176 _ 0,13681 

0,2440 0,3340 0,17577 

0,2886 0,34835 0,22124 

0,3362 0,3610 0.27373 

0,3877 0,3723 0,33373 

0,4424 0,3824 0,40174 

0,5005 0,3915 0,47830 

0,5619 0,3998 0,56379 

0,6265 0,4072 0,65882 

0,6948 0,4140 0,76386 

0,7662 0,4203 0,87940 

0,8411 0.4260 1,00594 

0,9192 0.4314 1,14399 

1,0008 0,4363 1,29404 

1,0856 0,4408 1,45660 

1,1738 0,44505 1,63440 

1,2654 0,4490 1,82120 

1,3603 0,4527 2.02343 

1,4585 0,4562 2,24182 


Если теплоемкости веществ, участвующих в реакции, вы¬ 
ражены с помощью интерполяционных формул 

с р = а-\- §Т + тГ 2 

или 

с. = а + р7’+Т , 7'~ 2 , 


248 л 

Т- А5 &8 - + ЛМР + Мф Т + лиД7'), 


до^ 

~т~ 
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где Да, Др, Ау и Ау' — алгебраические суммы коэффициентов 
в температурном ряду теплоемкости. 

Величины Ліо, Лі ь Лі 2 , Лі_ 2 для различных температур рас¬ 
считаны Темкиным и Шварцманом и сведены в табл. 3. 

Приведем пример расчета по методу Темкина и Шварцма¬ 
на. Используя данные предыдущего примера, рассчитаем 
ДОУдао для реакции (8). 

По Темкину и Шварцману 

д пО Д/Ч® 

ТЙЙГ “ -ляг - А5 2» - + м г др + М,М + ЛІ_ 2 Д Т '). 

Для Т = 1000 К 

7И 0 = 0,5088; М х =0,2463-10 3 ; Лі_ 2 = 0,2783-10 -5 . 

Для рассматриваемой реакции 

Да = -5,03; Др = 0,83-10~ 3 ; Д т ' = 1,71 -10 5 . 
Следовательно, 

ДШ 2 = ~ 41 . 46 - [°- 5088 (- 5 > 03 ) + 0,2463-0,83 + 

+ 0,2783 -1,71] =68 — 41,46 — (— 2,56 + 0,20 + 0,47) = 28,43; 
ДО'/ооо = 28 430 кал. 


Расчет по спектроскопическим данным 

Из постулата Планка при нуле градусов по шкале Кель¬ 
вина 5о°=0. Это означает, что при абсолютном нуле Со 0 — 
— #о°, следовательно, появляется возможность рассчитывать 
абсолютные значения термодинамических функций. 

В самом деле, \ 

0° = Щ-Т8° Г 
но 

н 0 = н° + ] Ср йт-, 

о 


Значит, 
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Таким образом, и энтальпия вещества, и его изобарный по¬ 
тенциал являются однозначными функциями теплоемкости. 

В последние 30 лет спектроскопические исследования ве¬ 
ществ в газообразном состоянии обеспечили возможность 
наиболее точного определения их теплоемкости и изменения 
этой величины с температурой. Это позволило рассчитать 
термодинамические функции и составить таблицы их значе¬ 
ний для многих веществ. 

С помощью методов статистической термодинамики, осно¬ 
вываясь на спектроскопических данных, можно непосред¬ 


ственно рассчитать теплоемкость вещества, его энтропию и 
приведенные значения термодинамических функций: приве¬ 
ла — #о° 

леннѵю энтальпию НІ — -~- и приведенное значение 


изобарного потенциала 



отсюда 

0° = #о° - ТФ* Г 

Соответственно 

АО° т = АН 0 °-ТАФ* т , (9) 

а так как 

д а° т 

\пК р = дг-г > 
то 

Д Ф* г Д# 0 ° 

1п К р = яТ • 

Величины АФт* и АН 0 ° для химической реакции подсчи¬ 
тываются по общему правилу: суммируются значения термо¬ 
динамической функции для продуктов и вычитается сумма 
значений функции для исходных веществ (с учетом стехиомет¬ 
рических коэффициентов). 

Тепловой эффект реакции при абсолютном нуле Л#о° 
определяется либо из спектроскопических данных, либо из 
термохимических данных для любой температуры с помощью 
уравнения Кирхгофа 

т 

АН о° = АН 0 - | АСрйТ. 
о 

Поскольку для большинства веществ известны энтальпии 
при 298 К, чаще всего исходят из теплового эффекта при 
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этой же температуре. Кроме того, температурный интеграл 
легко рассчитывается из значений Н т ° — #о°, приводимых 
в таблицах спектроскопических данных. Таким образам, 

Значения Н 0 ° вычислены для важнейших газов и приво¬ 
дятся вместе со значениями Ф т * в таблицах. Что касается 
конденсированных фаз, то для них расчет А Н 0 ° может вы¬ 
звать затруднения, которых можно избежать, если в расчет 
ввести вместо функции Ф т * функцию Фт*'*, причем 

/->0 г /0 

и? у* • 



Значения первого приведенного изобарного потенциала 
Фт* численно определяют статистическим методом непосред¬ 
ственно из молекулярных параметров, значение же второго 
приведенного потенциала вычисляют с помощью первого из 
соотношения 


ф** __ ф* 


н, 


298 


Яо« 


При использовании в расчете второго приведенного потен¬ 
циала отпадает необходимость в использовании ДЯ 0 ° и вме¬ 
сто него в расчет вводится тепловой эффект при 298 К Л// 2 ° 98 , 
который легко определяется из табличных данных. Расчет 
изменения изобарного потенциала и константы равновесия 
химической реакции ведется по формулам: 

до« = -гдф;* + д#« 98 ; 

При наличии табличных значений Ф Т * и Ф т ** для веществ, 
участвующих в реакции, расчет А С т ° таким способом очень 
удобен, ибо осуществляется простыми алгебраическими дей¬ 
ствиями (сложение или вычитание) и не требует применения 
интегрирования. К сожалению, таблицы термодинамических 
данных, рассчитанных из спектроскопических величин, не 
охватывают всех веществ и зачастую разбросаны в много¬ 
численных источниках, что затрудняет применение этого 
метода. 

Привёдем пример расчета по спектроскопическим данным 
для реакции (8). Расчет стандартного изменения изобарного 

2<? 



потенциала по формуле (9) требует знания первого приведен¬ 
ного потенциала для каждого участника реакции 

ЗіСІ™ 3 + 2Н™ 3 = 5і тв + 4НСГ аз . 

Кроме того, для расчета теплового эффекта реакции при 
абсолютном нуле А#о° необходимо знать тепловой эффект 
при 298 К Д// 2 ° 98 и значения функции Я 2 ° 98 —Я 0 ° для каж¬ 
дого участка реакции. Из справочника по* термодинамиче¬ 
ским свойствам, рассчитанным по спектроскопическим дан¬ 
ным, отыскиваем необходимые сведения: 


Вещество 

- * 

ф 1000’ 

кал/моль-град 

^298 — Д>°. 
кал/моль 

Н™ 3 

32,73 

2019 

5іС1™ 3 

86,32 

4633 

8і тв 

6,52 

770 

НСІ™ 3 

46,15 

2065 


Величину А // 2 °98 рассматрива емой реакци и заимствуем 
из предыдущего примера (Д# 2 ° 98 =^Ш) 0(Г ка л)А- Как и в пре¬ 
дыдущих примерах, будем вести расчет при 1000\. 

Учитывая стехиометрию реакции, производим расчет 
дф шо и д №- я о°): 

Аф* юо = 6,52 + 4 ■ 46,15 - 86,32-2-32,73= 191,12— 151,78 = 

= 39,34 кал/град; 

Д# в ° = 68000- [770 4-4-2065 -4633 -2-2019] =67 770 кал. 
Отсюда 

АО® 0(10 = — 39,34 -1000 + 67 770 = 28 430 кал. 

§ 5. Оценка термодинамической возможности 
протекания реакции 

Для типичных условий протекания химической реакции 
(7' = сопз1, Р = СОП5І) термодинамическая возможность про¬ 
текания реакции определяется знаком изменения изобарного 
потенциала реакции АО. Используя данные, приводимые 
в справочниках по термодинамическим свойствам веществ, 
можно рассчитать стандартное изменение изобарного потен¬ 
циала реакции АО 0 (для условий, в которых парциальное 
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давление любого газообразного участника реакции равно 
1 атм). Возникает вопрос, в какой мере стандартное измене-' 
ние изобарного потенциала реакции АО 0 мол&т служить кри¬ 
терием для определения термодинамической возможности 
протекания реакции? 

Анализ уравнения изотермы химической реакции (стр. 5) 
показывает, что знак АС 0 будет совпадать со знаком АО 
в двух случаях: 

1) парциальные давления газообразных участников реак¬ 
ции мало отличаются о*г давления в 1 атм; 

2) АО 0 велико по абсолютному значению (практически 
больше 10 ккал). 

В остальных случаях к оценке возможности протекания 
реакции по знаку АО 0 необходимо подходить с осторожностью 
и по возможности определять знак истинного изменения изо¬ 
барного потенциала. 

Пример. Будет ли окисляться серебро при температуре 
25° С в среде с парциальным давлением кислорода ІО -5 атм, 
если стандартное изменение изобарного потенциала реакции 

4А§: + 0 2 = 2А^О 

при этой температуре составляет 

АО° т — — 5,3 ккал? 

Согласно уравнению изотермы химической реакции 

ДО г = ДО° + ЯЛп —. 

Р о г 

Для наших условий 

Д(?298 = — 5300 — 4,575 ■ Т 1^ р 0> 

или 

ДС?298 == — 5300 — 4,575-298 (— 5). 

Следовательно, 

Д0 2Я8 = - 5300 + 6820 = + 1520 кал 
и окисление серебра термодинамически невозможно. 

§ 6. Расчет характеристик полноты протекания реакции 

Анализ любого химического процесса заключается нс 
только в установлении направления протекания реакции в за¬ 
данных условиях, но и в выявлении степени ее завершенно¬ 
сти. Термодинамической характеристикой, определяющей 
степень полноты реакции, глубину превращения исходных 
веществ в продукты, является константа равновесия химиче- 
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ской реакции. Знание этой характеристики позволяет рас¬ 
считать величину (степень превращения), непосредственно 
связанную с количеством вещества, вступившего в химиче¬ 
скую реакцию или образовавшегося в ходе ее протекания. 
В случае распада одного вещества эту величину часто назы¬ 
вают степенью диссоциации. 


Расчет степени диссоциации химического соединения 

Рассмотрим расчет степени диссоциации на примере реак¬ 
ции термического разложения тетрайодида кремния 

ЗУ™ 3 = 5і тв + 2^ аз . 

Определим степень диссоциации при температуре 1500 К и 
давлении 0,1 атм, если при этих условиях константа равнове¬ 
сия К р —1. Константа равновесия рассматриваемой реакции 

р\ 

Перейдем к числу молей 

„ 2 

К — р -Д— 

р п 8іі іп 

Если в реакцию вступает п' молей тетрайодида, а исход¬ 
ное число молей равно п 0 , то 

Л8и,=Ло-л'; 

2 п = я„ — п' + 2 п' — п 0 + п '. 

Тогда 

К —р -_ 

Р Г ІЛо-Я'Ж + П') 

Вынесем п 0 за скобку и обозначим степень диссоциации 
через 


Тогда 


кГ = Р 

Л Р г (1-«)(!+а) 


К —Р 4з ° ■. 

п Р 1 — а2 


Отсюда 




\р_ 

+ Кр 
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Для условий протекания нашей реакции 

“ = Кгахтт=°. 845 ; 

а = 84,5%. 

Расчет степени превращения 

Сначала рассмотрим расчет степени превращения для 
простого процесса, включающего только одну химическую 
реакцию. Определим степень превращения трибромида бора 
по реакции 

ВВг™ + 3/2 Н = В тв + ЗНВт газ . 

Условия проведения реакции: температура 1000 К, концентра¬ 
ция трибромида бора в смеси с водородом— 1 мол- %, дав¬ 
ление 1 атм. Константа равновесия реакции при температуре 
1000 К равна 0,01. 

Константа равновесия рассматриваемой реакции 


РвъцРщ 


Перейдем к числу молей 


^ВВГз-Ущ 


й ВВг 3 «Щ (Я")'' 2 ' 


Поскольку общее давление Р равно 1 атм, в дальнейшем его 
можно не учитывать. Если исходное число молей ВВг 3 равно 
«о.ь число молей ВВг 3 , вступивших в реакцию, — п', а исход¬ 
ное число молей водорода — п 0> 2 , то равновесные числа молей 
будут: 

П НВг~^ П > ^ВВГЗ^^О, 1 «> 

«На ~ «О, 2 3/2я ; 

ІЙ = з п! + п„ л — п' + л 0 2 -§- п'. 

Так как исходная газовая смесь является сильно разбавлен¬ 
ной (УѴ ВВгі =0,01), то 


Поэтому 


« 0 . 2 » П 0,1 И «'• 


«Н а ~ «0, 2* ^« ~ «0, 2- 
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Таким образом, 


К *7 (я ) 3 

р («о. і - «') («о, 2 ) 2 ' 

Вынесем Ло,і за скобку 

К 27 (я') 3 

Р «0,1 (! -Л'/Ло.і) л0.2 ’ 

Тогда, учитывая, что 

«7«о,і = 4, 

где т) — степень превращения ВВг 3 , получаем 

д- __ 27т , 3 _ 

Р (1 -Ч) («0,2/«0,і) 2 

или 

271 )<Вг а 

А и ” — ■« 1 • 

р 1 — Т, 

Отсюда получаем следующее уравнение для расчета степени 
превращения: 

27іѴ ВВг 

1 -4 = . щ гі 3 . 

Это уравнение проще всего решить графически, обозначив 


1 - ѵ) = у,; 


ввг., з 

яГ -71 = 


и определив пересечение этих кривых. 

Для условий проведения нашей реакции 

27Л г ввг 27 -(10 2 ) 2 

-Д ° Гз . = -._ і. — О 97 

К р 0,01 • 

Таким образом, степень превращения найдется из совмест¬ 
ного решения уравнений: 

1 -7] = у, ; 0,27т) 3 — у г - 

Из графика на рис. 6 находим г) = 84,3 %. 

Теперь рассмотрим пример расчета степени превращения 
для более сложного процесса, включающего несколько хими¬ 
ческих реакций. Процесс получения кремния восстановлением 
тетрахлорида кремния водородом характеризуется одновре¬ 
менным протеканием двух реакций: 

ЗіСІ™ 3 + 2Н г 2 аз = 5і тв + 4НСГ аз ; 

ЗІСІ™ + 5і тв = 23іС1 г 2 аз . 
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Для того чтобы найти степень превращения тетрахлорида 
кремния в кремний, необходимо записать константы равно¬ 
весия этих реакций, выразить их через степень превращения 
и решить совместно систему из двух уравнений. 

Определим степень превращения тетрахлорида кремния 
в кремний при температуре 1500 К, если константы^ равнове¬ 
сия указанных реакций соответственно равны: /Сі = 0,204, 
/С 2 = 0,135, давление—1 атм, а исходная концентрация тет¬ 



рахлорида в смеси с во¬ 
дородом — 2 мол - %. 

Запишем константы 
равновесия рассматривае¬ 
мых реакций: 

Л 

ѵ /’нсі . 

- 2~* у 

Р5ІС\,Рн 2 

Iгг ^5ІСі 2 

/*зІсй 


Перейдем от парци¬ 
альных давлений к числу 
молей: 


К\ — Р - 

«ЗіСі, «Н," л 


Рис. 6. Определение степени превраще¬ 
ния бромида бора в бор 


Кг 


= Р 


„2 

«5ІС1., , 

«ЗІСЦ^ ’ 


причем по условию Р— 1 атм. 

Если обозначить исходные числа молей тетрахлорида 
кремния и водорода через н 0 д и Яо.г, число молей 5іСЦ, 
вступающих в первую реакцию, через х, а во вторую через у, 
то равновесные числа молей будут: 

«на = 4х ; «зісі, “ «о,і ~~ х У > 

«8ІС1* — «н, “«0,2 

= п 01 — х — у + я 02 — 2х + 2у == ■ .: 

= «ОД + «0,2 — Зх + У . 


Так как по условию задачи исходная газовая смесь является 
разбавленной, то 


«0,2 > «о,н х и У- 

Поэтому 

«П а “ «0,2’ ^« ~ «0,2 
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и тогда 


_ 256** 

( п 0,\~ х -У) 4.2 ' 

К _ 4 У‘ 

2 («од -х-У) « 0,2 ’ 

Вынося /г 0 д за скобку и учитывая, что 

*/«0Д = Ѵ Й> 

где ці — степень превращения 5іСЦ в кремний; 

У/«од = 

где -Ц 2 — степень превращения ЗіСЦ в дихлорид кремния, 
получаем: 


К х = 


(1 — К)! — 1)2) («0,2/«0Д )'’ 

_ 4і)г 2 _ 

(1-1)1 - 1)2І «0,2/«0Д ' • 


Так как 


1 1 


«ОД ^БіС-І, 


К д 2561І.1 

1 1 — Т |1 — 1)2 ’ 


К 2 = ~.—^- N. 

г 1 — 1)1 - 1)2 


Решаем совместно полученную систему уравнений. Поделив 
первое уравнение на второе, выразим т )2 через тц 


__ 256гц* », 2 

К, 4Ѵ 


Отсюда 


Ч, = 8. 

Подставляем в первое уравнение 

К, (I - ъ - 8М -I), 2 ) = 256/ѴЧ 4 

и преобразовываем 

1 - ѣ == 256 V + 8УѴ 1)! 3 . 


3 
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Уравнение удобно решать графически, обозначив 

1 —т), = г,; 


|-Ѵ + 8Л^|/-|-ѵ = 2 2 . 


Для условия нашей задачи: 



О ,7 0,8 0,9 ІО 


266 -|- = 256 = 0,010; 

8Л^ ]/^- = 8 • 0,02 = 0,132. 

Графическое решение системы 
уравнений 

. 1 — Т)і = 2,; 

0,010т], 4 + 0.132Ѵ = 2, 

приводит (рис. 7) к значению гц — 
= 89 %. 

— і —— § 7. Расчет теплоемкости 


1 Для расчета температурной зави- 
Рис. 7. Определение сте- симости любой термодинамической 
пени превращения тетра- функции (теплового эффекта, изме- 
хлорида К Р®^“ ИЯ в к Р ем ‘ нения энтропии, изменения изобар¬ 
ного потенциала) необходимо иметь 
данные о теплоемкостях веществ, участвующих в реакции. 
Наиболее надежные данные по теплоемкости дают калори¬ 
метрические измерения. Справочные данные по теплоемкости 
задаются в виде температурных рядов. В последнее время 
наиболее широко используется уравнение, предложенное 
Майером и Келли, 

с р — »•+ Р-7’,+ ГТ 2 . 

Это уравнение хорошо передает температурную зависимость 
теплоемкости в широком интервале температур: при высоких 
температурах третье слагаемое очень мало и зависимость 
приобретает линейный характер, при низких температурах 
это слагаемое вносит существенный вклад в величину тепло¬ 
емкости и обеспечивает криволинейный характер зависимо¬ 
сти теплоемкости от температуры. 

Довольно часто используется уравнение типа 

с р = * + $Т-\- Л Т\ 
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Однако при высоких температурах оно приводит к большей 
кривизне температурной зависимости теплоемкости, чем это 
следует из экспериментальных данных. Кроме того, при 
использовании этого уравнения часто оказывается, что тре¬ 
тий член уравнения уТ 2 имеет отрицательный знак, т. е. при 
некоторой температуре теплоемкость проходит через макси¬ 
мум, что не имеет физического смысла. 

Если справочные данные по теплоемкости отсутствуют, 
то теплоемкости рассчитываются одним из следующих мето¬ 
дов. При наличии сведений о теплоемкости в ряду подобных 
веществ теплоемкости оцениваются по методу Лотье — Кара- 
петьянца (аналогично тому, как это было рассмотрено для 
теплоты образования и энтропии) или по методу химического 
подобия. 

Метод химического подобия 

Для расчета температурной зависимости теплоемкости по 
методу химического подобия пользуются температурной зави¬ 
симостью теплоемкости для вещества, химически подобного 
рассматриваемому. Примерами таких химически подобных 
веществ служат: 

СС1 4 и 8іС1 4 , ЗіС1 4 и 5іВг 4 , НС1 и НВг. 

Предположим, нужно найти теплоемкость вещества А 
в виде температурного ряда 

с г (.А)=‘ л +? А т+ит- 2 , ' 

и нам известен такой ряд для вещества В, химически подоб¬ 
ного веществу А, 

с р (В) = а а + М'+Т в 7"- 2 , 

сокращенно эти ряды можно записать в виде 

с р = 2X7"", ' 

где X — это а, р, у. Тогда 

с р {А)=Ы л Т”-, 
с р (В) — 2 Х В Г П . 

При некоторых температурах Т А и Т в : 

с р {А) — с р {В)\ 

Ср {А)-с р {В) = 0; 


3* 
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Вынося температуру Т А за скобку, получаем 

*ЩК-КЧІП)= о, 

откуда 

К-К[т в Іт А у. 

т 

Это означает, что аА = а в , так как п — 0; р А = (З в -=-,так 

* А 


как п= 1; ■у А = у в ("т^| > так как п — — 2. 

Как видно, для искомого вещества А надо знать хотя бы 
одно значение теплоемкости при какой-либо температуре. 
Тогда, зная с р (А) при Т А из известного температурного ряда 
для вещества В, находим Т в , при которой с р (В) = с р (А). 
Подсчитываем коэффициент подобия Т 0 <=Т В /Т А и опреде¬ 
ляем неизвестные коэффициенты: 

а А~ а в’ ^л == ?в'* г «’ ?л ” 2- 

В качестве примера рассмотрим определение теплоемкости 
газообразного ЗіСЦ методом химического подобия. Выби¬ 
раем подобное вещество, например ЗпСЦ. Теплоемкость газо¬ 
образного ЗпСЦ описывается следующим температурным 
рядом: 

с р = 25,57 + 0,20-10” 3 7' — 1,87- Ю 5 7'~ 2 . 

Известно, что теплоемкость газообразного ЗіСЦ при темпе¬ 
ратуре 300 К: 

с р (ЗіСЦ) = 21,67 кал/моль-град. 

Найдем температуру, при которой 

с р (ЗіСЦ) = с р (ЗпСЦ). 

Для этого в уравнение теплоемкости ЗпСЦ подставляем зна¬ 
чение 21, 67 и, решая уравнение, получаем 


, так как п ■■ 


Т 5п а- 218 К. 


Коэффициент подобия 


91 о 

“збо” ~ 0,7 12* 


Следовательно, 

“Цісі,” а 5 псі, = 25,57; 

Р5ю,= Рзпа/о = 0,20-0,712 = 0,14- ІО -3 ; 
Т5ІСЦ = Тзпсі, V ='87- (0,712)-* = 2,62 • 10 5 . 
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Таким образом, теплоемкость газообразного ЗіСЦ описы¬ 
вается следующим уравнением: 

с р = 25,57 + 0,14-ІО -3 /-- 2,62-10 3 Г“ 2 . (10) 

Как видно из приведенной ниже табл. 4, расхождение с экспе¬ 
риментальными значениями не превышает 5%. 

Таблица 4 


Расчет теплоемкости ЗіСЦ по методу химического подобия 


Температу¬ 
ра, К 

Теплоемкость, 

кал/моль-град 

Расхождение, % 

расчет по фор¬ 
муле (10) 

экспериментальн. 

значения 

300 

22,70 

21,67 

4,3 

500 

24,59 

24,02 

2,4 

1000 

25,45 

25,34 

0,4 


Теплоемкость жидкостей 

Теория теплоемкости жидкости разработана значительно 
хуже, чем теория теплоемкостей газов и твердых тел, и в на¬ 
стоящее время нет простых и точных количественных мето¬ 
дов определения этой величины для жидкого состояния веще¬ 
ства. При отсутствии экспериментальных данных теплоемко¬ 
сти жидкости оцениваются по методу атомных инкрементов, 
согласно которому 

где с р , і — атомная составляющая (инкремент) теплоемкости 
элемента; Пі — число атомов данного элемента в молекуле 
соединения. 

Атомные инкременты теплоемкости приводятся в справоч¬ 
ной литературе [4] или вычисляются из экспериментальных 
значений теплоемкости химических соединений. 

Пр и мер. Оценить теплоемкость жидкого тетрайодида 
олова ЗпСЦ, если известны теплоемкости^следующих жидких 
галогенидов: 

ЗіСЦ — 34,7 кал/моль • град; 

ЗпСЦ — 39,5 кал/моль-град; 

ЗП 4 — 35,1 кал/моль-град, 

а также атомный инкремент теплоемкости кремния, равный 
5,80 кал/г-атом-град. Поскольку 

с р (5і,Ц) = с р (8і) -}- 4 с р (Л), 
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і 


то 


Так как 
и 

то 


Ч(Л) = с р (51+) - с р (5і) = 35,1 — 5,8 = 
= 29,3 кал/г-атом-град. 

с р (ЗпС1 4 ) = с р (Зп) + Ас р (СА) 

с р (ЗіС1 4 ) = с р (Зі) + 4с р (СІ), 


с р (Зп) - с р $і) - 39,5 - 34,7 = 4,8, 

откуда 

^р(Зп) = (51) + 4,8 = 5,8 + 4,8 = 10,6 кал/г-атом-град. 


Таким образом, 

с р (ЗпЛ 4 ) == 10,6 4- 29,3 — 39,9 кал/моль - град. 


Экспериментальное значение теплоемкости жидкого 5п+ со¬ 
ставляет 40,1 кал/моль-град. 


Расчет теплоемкости твердых и газообразных веществ 
методами квантовой теории 

В классической статистической физике установлен закон 
о равномерном распределении энергии по степеням свободы, 
согласно которому на каждую степень свободы поступатель¬ 
ного или вращательного движения в среднем приходится 
кинетическая энергия, равная в расчете на одну частицу 
'ІъЬТ, в расчете на один моль вещества 1 / 2 РТ . Поскольку 
энергия гармонического осциллятора складывается из кине¬ 
тической и потенциальной энергии, причем в среднем они рав¬ 
ны между собой, то средняя энергия колеблющейся частицы 
в два раза больше, чем энергия частицы, совершающей посту¬ 
пательное или вращательное движение. Таким образом, на 
одну степень свободы колебательного движения по класси¬ 
ческой статистике приходится энергия, равная КТ. 

Закон о равномерном распределении энергии по степеням 
свободы позволяет легко рассчитать теплоемкость системы 
частиц, совершающих определенное движение. Рассмотрим 
определение теплоемкости идеального газа. Молекулы иде¬ 
ального газа участвуют в поступательном и вращательном 
движении. Поскольку теплоемкость, приходящаяся на одну 
степень свободы поступательного или вращательного движе¬ 
ния, равна; 

^ = тИ-5-« 7 ')„ = -5-*' 
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то теплоемкость идеального газа будет 

Су — 1/2КІ, 

где і — число степеней свободы молекулы. 

Таким образом, для одноатомного идеального газа 


і 




Для двухатомного идеального газа 


/=3+2=5" 



Для трех- и многоатомного идеального газа необходимо учи¬ 
тывать конфигурацию молекулы. Для линейных молекул 


+ ^'вращ —-3+2 



Для нелинейных молекул 

^ == ^поетуп + ^вращ == 3 + 3 = 6, 



Применение закона о равномерном распределении энер¬ 
гии к твердым кристаллическим веществам позволяет теоре¬ 
тически вывести закон Дюлонга и Пти, установленный экспе¬ 
риментально для теплоемкости металлов при комнатной тем¬ 
пературе. Рассматривая узлы кристаллической решетки твер¬ 
дого тела как трехмерные гармонические осцилляторы и учи¬ 
тывая, что энергия одномерного осциллятора в среднем 
равна кТ (в расчете на одну частицу) или КТ (в расчете на 
один грамм-атом вещества), получаем: один грамм-атом 
твердого вещества обладает энергией ЗКТ. Соответственно 
теплоемкость твердого тела 

с ѵ = ЪК = (з кал/моль-град (25 Дж/моль -град). 


Классическая теория теплоемкости хорошо описывает теп¬ 
лоемкости одноатомных реальных газов и многих простых 
кристаллических тел при комнатной температуре. В осталь¬ 
ных случаях она вступает в резкое противоречие с опытом, 
устанавливающим наличие температурной зависимости теп¬ 
лоемкости. Это обусловлено неприменимостью законов клас¬ 
сической физики к описанию явлений атомного масштаба, 
которые подчиняются особым, квантовым закономерностям. 


Квантовый закон распределения энергии 
по степеням свободы 

В основе квантовой теории лежит представление об испу¬ 
скании и поглощении энергии отдельными неделимыми пор¬ 
циями — квантами. В соответствии с этим энергетический 
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спектр каждого вида движения есть набор отдельных, 
дискретных энергетических уровней. Переход микрочастицы 
с одного уровня на другой происходит при поглощении или- 
испускании кванта энергии. В случае поступательного или 
вращательного движения расстояние между отдельными 
энергетическими уровнями настолько мало, что спектр этих 
видов движения для не слишком низких температур можно 
считать практически непрерывным. Следовательно, кванто¬ 
вое рассмотрение поступательного и вращательного движе¬ 
ний совпадает с классическим и для этих движений сохра¬ 
няется закон о равномерном распределении энергии по сте¬ 
пеням свободы. Расхождение между классической и кванто¬ 
вой теориями наблюдается только при рассмотрении колеба¬ 
тельного движения, для которого закон о равномерном рас¬ 
пределении теряет свою силу. 

Для того чтобы вывести квантовый закон распределения 
энергии по колебательным степеням свободы, рассмотрим 
систему, состоящую из N осцилляторов, среди которых рас¬ 
пределено УѴ 0 квантов, каждый из которых несет энергию йѵ. 
Полная энергия нашей системы 

Д=ЛУЬ. (И) 

Энергия по осцилляторам распределена так, что одни из 
них обладают энергией йѵ, другие — энергией 2 йѵ, третьи — 
энергией Зйѵ и т. д. Наша задача — найти энергию, прихо¬ 
дящуюся на частоту ѵ, т. е. энергию осциллятора, колеблю¬ 
щегося с частотой ѵ. 

В равновесном состоянии распределение квантов по осцил¬ 
ляторам должно удовлетворять формуле Больцмана 

5 = к \пѴУ, (12) 

где ѴУ — термодинамическая вероятность системы, равная 
числу возможных микросостояний системы. В нашем случае 
микросостояние системы —- это определенное распределение 
квантов по осцилляторам. Общее число возможных микросо¬ 
стояний равно числу способов, которыми УѴ 0 квантов можно 
разместить по N осцилляторам. 

Задача о нахождении числа способов, которыми УѴ 0 кван¬ 
тов можно разместить среди N осцилляторов, эквивалентна 
задаче размещения УѴ 0 точек по N клеткам, разделенным пе¬ 
регородками. Число возможных микросостояний равно числу 
способов размещения точек по ячейкам, причем микросостоя¬ 
ние меняется, если меняется число точек в какой-либо ячейке. 
Каким образом можно изменить количество точек в ячейке? 
Это можно осуществить с помощью перестановки (обмена 
местами) любых двух элементов системы. В нашей системе, 
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состоящей из точек и перегородок (рис. 8), возможны сле¬ 
дующие перестановки: 

а) перестановка точек; 

б) перестановка перегородок; 

в) перестановка точки и перегородки. 

Поскольку перегородок на единицу меньше, чем клеток, 
то общее число перестановок всех элементов (и точек, и пе¬ 
регородок) равно (Ы 0 + N — 1)! Но перестановка самих то¬ 
чек или перегородок не изменяет микросостояния системы, 
поэтому их надо исключить из общего числа перестановок. 
Поскольку число перестановок точек равно Л/ 0 !, а число пере¬ 
становок перегородок равно (—— 1)!, то число перестановок, 



Рис. 8. Подсчет числа микросостояний системы из 
Л ; о квантов, распределенных между N осцилля¬ 
торами 

Таким образом, число микросостояний системы из А г 0 
квантов и N осцилляторов 

(іѴо-НѴ-1)! _.(А *о + Щ 

™ ~ Щ(Ы~ 1)1 лум ’ 

так как 

К и УѴ» 1. 

Подставляя полученное выражение в формулу для энтропии 
(12) и используя формулу Стирлинга 

1пх! =л;1п(х — 1), 

которая для больших х трансформируется в 

1п х\ =ЛІПЛ, 

получаем 

~ = (УѴ 0 + УѴ) 1п (К + УѴ) - УѴ 0 1п УѴ 0 - УѴ 1п N. 

Из объединенного уравнения первого и второго начал термо¬ 
динамики 

таз=аЕ+ра ѵ 
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следует, что 


/ \ л 

[дЕ ) у~ Т 


Энергия нашей системы дается уравнением (11), подста¬ 
новка которого в последнюю формулу приводит к выра¬ 
жению 

1 / д8 \ _ 1 

Лѵ \дМ 0 )ѵ Т 

Но 

ж)- 


Следовательно, 


откуда 


*Ч'+т)~г. 

1п ( 1 +ж) в ‘ё г 


1 I Л _ р ьт 

N. ~ е * 


Последняя формула показывает, сколько квантов энергии 
приходится на один осциллятор. Поскольку каждый квант 
несет энергию Лѵ, то энергия осциллятора с частотой ѵ 


= ^-Ь = 
N п 


кТ 

е — 1 


Полученное выражение является математическим выраже¬ 
нием квантового закона распределения энергии по колеба¬ 
тельным степеням свободы. 


Теория теплоемкости Эйнштейна 

С появлением квантовых представлений Эйнштейн (1907 г.) 
сразу же предложил рассматривать твердое тело как сово¬ 
купность трехмерных осцилляторов, колеблющихся с одной 
и той же частотой ѵ. Используя выражение для энергии 
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одного осциллятора (13), для теплоемкости одной колеба¬ 
тельной степени свободы, получаем 


ѵ -(■$•). 


Лѵ 

Тт Лѵ 1 
е к * Г 2 

см ■ 


Умножим и разделим это выражение на постоянную Больц¬ 
мана 

Лѵ 

( Лѵ \2 кТ 

. \Тт) е 

с,, — к , —г=-. 


Лѵ \3 • 

лг 


Учитывая, что теплоемкость трехмерного осциллятора будет 
в три раза больше, и ведя расчет на один грамм-атом твер¬ 
дого тела (число осцилляторов равно числу Авогадро), по¬ 
лучим 


ЗЛ \к 


Лѵ 

Лѵ \ 2 Тт 

ит ) 2 


(-&) 


Так как = то 


с ѵ = ЪЯ. 


\" 

-і) 

(*) 

Лѵ 

2 кТ 

1 е 

/ Лѵ 

\ 2 

\е кТ - 



Если обозначить Нѵ/кТ — х, то последнюю формулу можно 
представить в более компактном виде: 

«Ѵ-ЗК-р^. (15) 

Уравнения (14) и (15) называются уравнениями Эйнштейна. 

Проведем анализ уравнения Эйнштейна для разных тем¬ 
пературных областей. Рассмотрим область высоких темпера¬ 
тур, в которой значение аргумента х мало. Тогда, разлагая е х 
в ряд по степеням х и ограничиваясь вследствие малости х 
двумя первыми членами 

\+х, 

получим 

с _ 3 Р-^1.Нг*.)- 
4 V ОА (1+л:-1) 2 
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или 


с ѵ — 3/?(1 +*)• 

Нол: — малая величина, следовательно, с ѵ —Зі?. Таким обра¬ 
зом, в области высоких температур формула Эйнштейна 
переходит в классическую формулу для теплоемкости твер¬ 
дого тела. 

Рассмотрим теперь область низких температур: Т-*■ 0; 
х-уоо. Тогда в формуле 

г — яр хіе * 

е х ->-оо и по сравнению с ней можно пренебречь единицей 
в знаменателе. Получаем 

^ = з4г = 3«Л- 

Экспоненциальный множитель оказывается доминирующим, 
вследствие чего при низких температурах с ѵ убывает с тем¬ 
пературой по экспоненте 

Н'* 

с^е~ х = е~ ш . 

Раскрыв неопределенность по правилу Лопиталя 
Пт ЪКх*е- х = 3/? Ііш — 3/? Нт -Д-, 

Х-юо ДГ-*-оо в ЛГч-оо ^ 


получим 


Нт с„ = 0, 


что хорошо совпадает с экспериментальными данными по 
теплоемкости при низких температурах. 

Величина Ѳ = кѵ/кТ носит название характеристической 
температуры. В самом деле ее размерность 

гді — ММ эрг-с-1/с _ . 

I ‘ [к] эрг/град Р 

Вводя характеристическую температуру в уравнение Эйнш¬ 
тейна, получаем 


с = -. 

К ( е 9/Г_1) 2 

Отсюда видно, что теплоемкость есть однозначная функция 
переменной Ѳ/Г: 

Эта функция носит название функции Планка — Эйнштейна. 
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Согласно модели Эйнштейна каждый атом вещества рас¬ 
сматривается как гармонический осциллятор, совершающий 
колебания независимо от колебаний других атомов. Если го¬ 
ворить о кристалле, то его спектр колебаний атомов согласно 
модели Эйнштейна имеет вид, представленный на рис. 9. 
Здесь по оси абсцисс отложена частота колебаний, а по оси 
ординат — функция распределения колебаний по частотам, 
т. е. число колебаний, приходящихся на интервал частот от ѵ 
до ѵ + кѵ. Такое допущение более справедливо для газов, 
нежели для твердых тел, поскольку взаимодействие атомов 

Су 


> і 

Рис. 9. Спектр колебаний по Эйн- Рис. 10. Температурная за- 
штейну висимость теплоемкости кри¬ 

сталла: 

1 — теория Эйнштейна: 2 — эксперимент 



кристалла приводит к тому, что атомы колеблются с самыми 
разнообразными частотами. Поэтому теория теплоемкости 
Эйнштейна дает лучшие результаты в применении к газам. 

В случае кристаллов эта теория хорошо описывает факт 
резкого уменьшения теплоемкости при низких температурах, 
но ход этого спада оказывается отличным от эксперименталь¬ 
ного (рис. 10). По теории Эйнштейна спад имеет экспонен¬ 
циальный характер, опыт же дает зависимость, пропорцио¬ 
нальную Г 3 . Что касается высоких температур, то здесь тео¬ 
рия Эйнштейна дает хорошее совпадение с опытом, и 
расчет теплоемкости твердых тел при температурах, от ком¬ 
натной и выше, можно вести по этой теории. 

Расчет теплоемкости газов 


Энергия газовой молекулы складывается из энергий по¬ 
ступательного, вращательного и колебательного движений 
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Соответственно теплоемкость газа есть сумма теплоемкостей 
указанных движений 

с ѵ = 6’™ СТ У П -{- С В Р* Ш + с к °- кЪ . 

Если рассчитать спектр энергии поступательного движения 
по уравнениям квантовой механики, то окажется, что рас¬ 
стояние между двумя соседними уровнями настолько мало, 
что им можно пренебречь. Следовательно,, и классическая, и 
квантовая физика в случае поступательного движения при¬ 
водят к одному и тому же результату: спектр энергии посту¬ 
пательного движения непрерывен, и справедлив закон равно¬ 
мерного распределения энерТии по степеням свободы. Сле¬ 
довательно, для любой молекулы 

Поступательное движение и поступательная теплоемкость 
оказываются единственными, где классическая и квантовая 
теории совпадают. В случае вращательного движения кван- і 
тованием энергетического спектра пренебрегать уже нельзя. 
Соответственно величина вращательной теплоемкости будет ' 

зависеть от температуры. Однако расстояние между враща- ; 

тельными уровнями невелико и уже при сравнительно невы¬ 
соких температурах (несколько десятков градусов по Кель¬ 
вину) они начнут участвовать в поглощении энергии. Поэтому 
для не слишком низких температур, при которых энергия 
теплового движения значительно больше расстояния между 
двумя соседними вращательными уровнями, спектр энергии 
вращательного движения можно считать непрерывным, и для 
расчета вращательной теплоемкости применяют простую 
классическую теорию 

Следовательно, для линейных молекул 

а для нелинейных . 

серии = -!_/?. 

Таким образом, для практических расчетов теплоемкости, 
которые ведутся в интервале температур от комнатной и 
выше, температурной зависимостью вращательной теплоем- : 
кости можно пренебречь. Следует, однако, иметь в виду, что 
при очень низких температурах квантование вращательной 

46 


энергии приводит к уменьшению величины вращательной 
теплоемкости по мере уменьшения температуры, происходит 
как бы «замораживание» вращательных степеней свободы, 
выключение их из процесса поглощения тепла. 

Расчет колебательной составляющей теплоемкости ведется 
по уравнению Эйнштейна 

Лѵ 



где с е (Ѳ/Г) — функция Планка — Эйнштейна; т —- число 
колебательных степеней свободы. 

Значения функции Планка — Эйнштейна приводятся в 
табл. 5. 

Т аблица 5 


Значения функции Планка — Эйнштейна 


ѳ 

Т 

кал 

Ѳ 

Т 

<Ч> 

кал 

ц 

рви 

моль - град 

моль-град 

моль-град 

0 

1,986 

1,45 

1,673 

3,70 

0,707 

0,10 

1,983 

1,50 

1,659 

3,80 

0.672 

0,15 

1,981 

1,55 

1,633 

3,90 

0,637 

020 

1,979 

1,60 

1,612 

4,00 

0,604 

0,25 

1,976 

1,65 

1,592 

4,20 

0,542 

0,30 

1,974 

1,70 

1,570 

4,40 

0,484 

0,35 

1,967 

1,75 

1,549 

4,60 

0,431 

0,40 

1,960 

1,80 

1,527 

4,80 

0,383 

0,45 

1,952 

1,85 

1,505 

5,00 

0,339 

0,50 

1,945 

1,90 

1,483 

5,20 

0,300 

0,55 

1,938 

1,95 

1,461 

5,40 

0,262 

0,60 

1,928 

2,00 

1,439 

5,60 

0,232 

0,65 

1,918 

2,10 

1,393 

5,80 

0,204 

0,70 

1,908 

2,20 

1,348 

6,00 

0,178 

0,75 

1,896 

2,30 

1,302 

6,40 

0,136 

0,80 

1,884 

2,40 

1,256 

6.80 

0,103 

0,85 

1,871 

2,50 

1,210 

7,20 

0,077 

0,90 

1,858 

2,60 

1,164 

7,60 

0,057 

0,95 

1,844 

2,70 

1,119 

8,00 

0,0427 

1,00 

1,829 

2,80 

1,074 

8,40 

0,032 

1,05 

1,814 

2,90 

1,030 

8,80 

0,023 

1,10 

1,798 

3,00 

0,986 

9,20 

0,017 

1,15 

1,782 

3,10 

0,943 

9.60 

0,012 

1,20 

1,765 

3,20 

0,901 

10,00 

0,009 

1,25 

1,747 

3,30 

0,860 

И 

0,004 

1,30 

1,729 

3,40 

0,820 

12 

0,0017 

1,35 

1,711 

3,50 

0,781 

13 

0,0007 

1 40 

1,692 

3,60 

0,744 
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Число колебательных степеней свободы можно опреде¬ 
лить следующим образом. Если молекула состоит из я ато¬ 
мов, то она имеет Зя степеней свободы (каждый атом обла¬ 
дает тремя степенями свободы и все они сохраняются в мо¬ 
лекуле) . Из общего числа Зя степеней свободы на долю по¬ 
ступательного движения приходится 3 степени свободы, а на 
долю вращательного движения: 2 — для линейных молекул 
и 3 — для нелинейных молекул. Следовательно, число коле¬ 
бательных степеней свободы 

Зя — 5 — линейные молекулы; 

Зя — б — нелинейные молекулы. 

В общем случае каждой колебательной степени свободы 
соответствует своя собственная частота колебаний и, следо¬ 
вательно, своя характеристиче¬ 
ская температура. Однако в не¬ 
которых молекулах могут совер¬ 
шаться несколько колебаний с 
одной и той же частотой. В та¬ 
ком случае говорят о вырож- 
01 дении колебаний, а число ко¬ 
лебаний с одинаковой часто¬ 
той называют кратностью вы¬ 
рождения. Вырождение колеба¬ 
ний характерно для симметрич- 
Рис. 11. Молекула тетрахло- ных молекул. В этом случае 
рида кремния слагаемые в выражении для ко¬ 

лебательной теплоемкости (16) 
входят с множителем, численно равным кратности вырож¬ 
дения. 

Произведем подсчет числа колебательных степеней сво¬ 
боды на примере молекулы газообразного четыреххлористого 
кремния, имеющей форму тетраэдра (рис. 11). Общее число 
степеней свободы этой молекулы 

і — Зя = 3-5 — 15. 

Число колебательных степеней свободы 

т — Зя — 6 = 9. 

Оказывается, из этих 9 колебаний одно невырожденное, одно 
двукратно вырожденное и два трехкратно вырожденных. 

Частоты колебаний молекул находятся из инфракрасных 
спектров поглощения и спектров комбинационного рассеяния 
молекул. Данные о частотах колебаний можно найти в спра¬ 
вочной литературе. 
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Следует отметить, что в справочниках приводятся не ча¬ 
стоты колебаний, а волновые числа (хотя последние име¬ 
нуются частотами колебаний). Волновое число есть обратная 
величина длины волны ѵ— 1Д. Учитывая связь между дли¬ 
ной волны, частотой колебаний и скоростью их распростра¬ 
нения (скорость света), получаем: 



Принимая для скорости света значение 3• ІО 10 см/с, получаем 
следующую формулу для расчета характеристической темпе¬ 
ратуры колебания с волновым числом ѵ: 

Ѳ = 1,413ѵ, 

если волновое число выражено в обратных сантиметрах. 
Переход от с ѵ к с р для газов осуществляется по уравнению 
Майера 

с р = с ѵ -(- К. 

Пример. Рассчитать теплоемкость газообразного тет¬ 
рахлорида кремния ЗіСЦ при 1000 К по теории Эйнштейна. 

Суммарная изохорная теплоемкость газа складывается 
из теплоемкостей поступательного, вращательного и колеба¬ 
тельного движений 

С __ ^.поступ ^вращ _|_ ^колеб. 

причем 

'СРОСТУ" = - 3- Я; С"Р аЩ = 4 Я 

(молекула ЗіСЦ относится к типу нелинейных молекул и 
обладает тремя степенями свободы вращательного дви¬ 
жения) . 

По теории Эйнштейна 

с Г еб = (-?-)’ 

1 

где 6 — число колебательных степеней свободы молекулы 

31СІ4. 

Из справочной литературы [6] находим характеристики 
колебательного движения молекулы ЗіСЦ — волновые чис¬ 
ла, см -1 : 

^ = 424; Ѵ 2 = 150(2); 7, = 621(3), 7 4 = 221 (3). 

Цифры в скобках указывают кратность вырождения. 
Таким образом, колебательное движение молекулы ЗіСЦ 
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описывается четырьмя характеристическими температурами. 
Соответственно ; 

С леб = (т-) + 2с г (4-) + Зс, (-^) + Зс е (А-) • 

Рассчитываем частоты колебаний по формуле 
ѵ = с7=3-10 10 7-Ъ-. 

С 

Получаем 

ѵ, = 1,272-ІО 13 ; ѵ 2 = 0,300-ІО 13 ; ѵ 3 = 1,863-ІО 13 ; 
ѵ 4 = 0,663- ІО 13 . 

По формуле 

О = 4,71 •10“" ѵ, 

где 4,7-10 ~ и = И/к, определяем характеристические темпе¬ 
ратуры, К: 

Ѳ, = 598; Ѳ 2 = 141; Ѳ 3 = 877; Ѳ 4 = 312. 

Рассчитываем аргументы функции Эйнштейна и по табл. 5 
определяем соответствующие значения колебательной тепло¬ 
емкости: 

ЩО = 0,598 = 1,927; 

1^0 = 0,141 Се =1,980; 
тш = 0.877 Се =1,864; 


0,312 Се =1,972. 


Таким образом, 


скол'б = , >927 + 2 . 1,980 + 3-1,864 + 3-1,972 =17,39 . 

Суммарная изохорная теплоемкость 5іСЦ 

с, = 3 + 3+17,39 = 23,3 9і ^^. 

Изобарная теплоемкость 

с р = с, + /?^26,39 + 2 = 25,39-^^. 

Экспериментальное значение изобарной теплоемкости ЗІСІ4 
при 1000 К составляет 25,62 кал/моль • град. 


Теория теплоемкости Дебая 

Физически необоснованное предположение о наличии 
одной частоты колебаний у всех атомов твердого тела и не¬ 
согласии в характере поведения теплоемкости вблизи абсо¬ 
лютного нуля привели к возникновению новой теории тепло¬ 
емкости твердых тел — теории Дебая (1919 г.). По Дебаю 
различные атомы колеблются в твердом теле с разными ча¬ 
стотами и из-за большого числа атомов спектр собственных 
частот можно считать практически непрерывным. Определе¬ 
ние теплоемкости в этом случае сводится к нахождению энер¬ 
гии системы атомов, колеблющихся с разными частотами. 

Согласно квантовой теории энергия, приходящаяся на 
одно колебание с частотой ѵ, равна: 


Полная энергия кристалла 


V * ЬТ 


где суммирование ведется по всем частотам колебаний кри¬ 
сталлической решетки. Поскольку по предположению Дебая 
спектр частот твердого тела считается непрерывным, содер¬ 
жащим все частоты от нулевой до некоторой максимальной 
частоты х т , суммирование можно заменить интегрированием 

ѵ т 

^=Г 

кТ 

ч) е — 1 

и 

где йг — число колебаний в интервале частот от ѵ до ѵ +]йх. 
Если рассмотреть функцию распределения колебаний по ча¬ 
стотам {(х) (число колебаний в единичном интервале частот 
от ѵ до ѵ + йх), то 

йг = / (ѵ) Ф*. 

Таким образом, расчет энергии кристалла сводится к нахож¬ 
дению функции распределения колебаний по частотам, дру¬ 
гими словами, к определению спектра колебаний кристалли¬ 
ческой решетки, для чего рассмотрим колебания его атомов 
в пространстве волновых векторов (все расчеты ведем на 
единицу объема кристалла). Как известно, 

1 _ 

X • 


ѵ = 
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В пространстве волновых векторов каждое колебание с ча¬ 
стотой ѵ изображается отдельной точкой. Число колебаний 
в интервале частот от ѵ до ѵ + (іѵ в пространстве волновых 
векторов соответствует числу состояний в интервале значе¬ 
ний волнового вектора от ѵ до ѵ + йѵ и будет равно объему 
шарового слоя (рис. 12) 

л йг = 4 лѵ 2 ^ѵ. 


Поскольку 



йч = — січ, 

V ' 


где ѵ — скорость распро¬ 
странения колебаний, то 
число колебаний в интер¬ 
вале частот от ѵ до ѵ + йѵ 
равно: 


Рис. 12. Подсчет числа колебаний 
кристаллической решетки твердого 
тела 


Так как каждый атом имеет три колебательных степени сво¬ 
боды, то полученное число надо утроить 


Если кристаллическая решетка содержит Л^ 0 атомов, то число 
колебательных степеней свободы равно ЗА/ 0 — б =* ЗА^о, по¬ 
скольку Л^ 0 очень велико (для одного грамм-атома ІѴ 0 ~ 
^6 -ІО 23 ). Следовательно, полное число колебаний 


} <йг= 5 /(ѵ)й?ѵ = ЗУѴ 0 


т . 

|-і^3ѵ^ = злѵ 


. Интегрируя, получаем 


52 


откуда 


Подставляем это выражение в уравнение (17): 

А-«***,, (18) 

откуда функция распределения колебаний по частотам 

X I \ о о 

/(^) = т^ ѵ 

*т . { 

выражается параболической зависимостью. Графически 

спектр колебаний кристаллической решетки по Дебаю имеет 
вид, изображенный на рис. 13; главный вклад в энер¬ 
гию вносят высокие ча- ,, , 
стоты, т. е. энергичные Тлѵ; 
кванты, а они возбуждают- і 

ся только при высоких тем- /| 

пературах. /1 

Подставляя формулу для / і 

числа колебаний в единич- / ; 

ном интервале частот (18) / і 

в выражение для энергии / | 

кристалла, будем иметь ! 


= -- I -- и-’. 

ч т 1 — Рис. 13. Спектр колебании кристалла 

Д е кТ _ і по Дебаю 

Введем новую переменную х = кѵІкТ, а также обозначим 
Ѳд=/гѵ то /&. Величина Ѳ с называется дебаевской температу¬ 
рой. Тогда выражение для энергии кристалла (в расчете на 
1 г-атом М 0 й = Д) будет 

.''я. г 

9 КТ А * 3 

4 ' О 

Часто его записывают в сокращенном виде 

Е = ЪЯТЭ [ ІвЛ , 


е х — 1 


называют функцией Дебая. 
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Интеграл, входящий в функцию Дебая, аналитически вы¬ 
числить нельзя. Исследуем его значение для двух предель¬ 
ных случаев. 

Рассмотрим высокие температуры, когда аргумент х— 
малая величина. В этом случае в формуле для функции 
Дебая (19) е х можно разложить в ряд по степеням х и 
ограничиться первыми двумя членами разложения: 

е* =* 1 + х. 

Тогда знаменатель обращается в х и функция Дебая преоб¬ 
разуется к виду 





хз Ь 0 /Т 
3 о 


следовательно, 

Е = с„ = ЗЛ, 


т. е. совпадает с классическим приближением. 

Рассмотрим теперь низкие температуры, когда аргу¬ 
мент х — большая величина. Разобьем участок интегрирова¬ 
ния от 0 до Ѳ 0 /Т на два интервала: от 0 до оо и от оо до Ѳв/7 



Первый интеграл 


хЫх п* 
е* - 1 “ 15“ • 

о 

Во втором интеграле, поскольку х велико, е х тоже велико 
и можно пренебречь в знаменателе единицей 

тг Г хЧ~*йх. 


При возрастании х сомножитель е~ х уменьшается гораздо 
быстрее, чем возрастает х 3 , поэтому при достаточно боль¬ 
шом х подынтегральное выражение становится пренебрежимо 
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малым и этим интегралом вообще можно пренебречь. Итак, 
при низких температурах 

Р_ЧР т 3 т -' 

п — эъі в з 15 


р _ 74 

15в 0 • 


При этом теплоемкость 

Эл 4 /? • 4 
с ѵ~ і5Ѳп 


Г 3 = 464,5 


пропорциональна третьей с ѵ 
степени абсолютной темпе¬ 
ратуры (закон 7-кубов Де¬ 
бая). 

Сравнивая кривую теп¬ 
лоемкости, рассчитанную по 
теории Дебая, с опытной 
кривой, следует отметить 
очень хорошее совпадение 
их в области очень низ¬ 
ких температур, при кото¬ 
рых теория Эйнштейна да¬ 
вала существенное расхож¬ 
дение с экспериментом. 

Так же, как и теория Рис. 14. Температурная зависимость 
Эйнштейна, теория теплоем- теплоемкости кристалла: 

КОСТИ Дебая хорошо совпа- /- теория Дебая; 2— эксперимент 

дает с опытом при высо¬ 
ких температурах. Небольшое расхождение теории Дебая 
с экспериментом наблюдается в области средних температур, 
примерно от 0,1 Ѳ с до 0,5 Ѳ п (рис. 14). Из опыта также сле¬ 
дует, что дебаевская температура Ѳв меняется с температу¬ 
рой, однако это изменение обычно не превышает 10%. 

Наличие температурной зависимости характеристической 
температуры' и небольшое отклонение от экспериментальной 
температурной зависимости теплоемкости свидетельствуют 
о приближенном характере теории Дебая. Это в первую оче¬ 
редь обусловлено отличием действительного спектра колеба¬ 
ний кристаллической решетки от параболического, дебаев¬ 
ского (рис. 15). • 

Однако, поскольку для определения теплоемкости произ¬ 
водится интегрирование функции распределения по частотам, 
это отличие сказывается очень мало на величине конечного 
результата. 
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Таким образом, в целом теория Дебая хорошо описывает 
температурную зависимость теплоемкости твердых тел. Необ¬ 
ходимо отметить, что Дебай развил свою теорию для изо¬ 
тропных тел, поэтому она хорошо описывает теплоемкость 
кристаллов с кубической симметрией, но несправедлива для 
анизотропных кристаллов (цепочечных или слоистых), на¬ 
пример для серы или графита. 

Главное достижение теории Дебая состоит в том, что она 
правильно описывает ход кривой теплоемкости при низких 
температурах. Это абсолютно точно выполняется для кри- 



1 - эксперимент; 2 -теория Дебая теристической температуры твер¬ 

дого тела 


сталлов изоляторов и полупроводников. У металлов расхож¬ 
дение теории Дебая с опытом наблюдается еще и при тем¬ 
пературах, меньших одного градуса по Кельвину, где стано¬ 
вится существенным вклад электронной составляющей тепло¬ 
емкости. 

Следует отметить, что закон Дюлонга и Пти, являющийся 
предельным случаем и теории Эйнштейна, и теории Дебая, 
лучше всего выполняется для типичных металлов. Например, 
при комнатной температуре теплоемкость меди равна 5,9, се¬ 
ребра— 6,0 и золота — 6,1 кал/г-атом. Закон Дюлонга и 
Пти не соблюдается для кристаллов, имеющих существенную 
составляющую ковалентной связи. Так, в четвертой группе 
периодической системы по мере перехода химической связи 
от ковалентной к металлической теплоемкость при комнат¬ 
ной температуре меняется так: алмаз—1,5; 5і — 4,7; Ое— 
5,3; 5п — 6,4 кал/г-атом-град. Мерой отклонения теплоем¬ 
кости простого вещества от закона Дюлонга и Пти является 
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характеристическая температура. Чем меньше величина 0 П , 
тем быстрее (с ростом температуры) теплоемкость достигает 
классического значения 3/?. Соответственно высокое значе- 
' ние характеристической температуры указывает на сущест¬ 
венное отклонение величины комнатной теплоемкости от клас¬ 
сического предела 3/? (рис. 16). Характеристические темпе¬ 
ратуры некоторых элементов приводятся в табл. 6. 

Таблица 6 

Температуры плавления и характеристические температуры 


некоторых веществ в кристалллическом состоянии 


Элементы 

Тпл. К 

е 0 , к 

Элементы 

Тпл. К 

е в , К 

А? 

1234 

221 

Мо 

2888 

377 

А1 

933 

390 

м 2 

62,7 

68 

А§ 

І090 

275 

N3 

370,7 

155 

А и 

1336,2 

178 

N6 

2741 

260 

В 

2498 

1362 

N6 

1297 

148 

Ва 

998 

116 

N1 

1726 

345 

Ве 

1557 

1031 

Ир 

910 

163 

Ві 

544.5 

116 

0 2 

54,7 

89 

Вг 2 

265,8 

86 

Оз 

3300 

400 

С (алмаз) 

4100 

1874 

Р (белый) 

317 

576 

Са 

1112 

230 

РЬ 

600,6 

87 

Се 

1070 

138 

Рб 

1825 

275 

Сй 

594,2 

221 

Рі 

2042 

225 

Сі 2 

171,9 

115 

Рг 

1208 

138 

Со 

1765 

386 

Рн 

913 

176 

Сг 

2148 

424 

кь 

311,8 

59 

Св 

301,8 

43 

Ке 

3433 

275 

Си 

1356 

310 

Р6 

2233 

350 

г 2 

50 

100 

К и 

2553 

415 

Ре 

1808 

373 

5с 

1812 

476 

Оа 

302,8 

240 

8і 

1700 

689 

Ое 

1209 

403 

8г 

1045 

148 

ОЙ 

1585 

155 

8Ь 

903,6 

200 

Н 2 

13,96 

91 

5п (белое) 

506 

170 

Ні 

2495 

213 

Та 

3271 

225 

н§ 

234,3 

92 

ТО 

2024 

100 

Но 

1734 

161 

Ті 

1941 

380 

І2 

386,6 

76 

Т1 

576 

96 

ІП 

429,8 

129 

11 

1404 

300 

Іг 

2716 

228 

V 

2178 

390 

К 

336,6 

100 

\У 

3635 

312 

Ьа 

1193 

135 

У 

1775 

214 

и 

454 

448 

2п 

692,6 

237 

Щ 

923 

330 

2т 

2123 

250 

Мп 

1517 

363 





Расчет теплоемкости твердых тел 

Теплоемкость простых твердых кристаллических тел рас¬ 
считывается по теории Дебая. Для этого определяют значе¬ 
ние аргумента Ѳ в /Г и по таблице дебаевской теплоемкости 

57 




находят соответствующее значение с ѵ . Теплоемкость слож¬ 
ных кристаллических тел определяют по формуле Неймана — 
Коппа, согласно которой теплоемкость химического соедине¬ 
ния равна сумме атомных теплоемкостей элементов 

і 

В кристаллической решетке соединения атом колеблется ' 
с частотой ѵ', отличной'от частоты колебаний в кристалли¬ 
ческой решетке простого вещества ѵ. Соответственно в раз¬ 
ных соединениях атом данного элемента имеет разные часто¬ 
ты. Поэтому в случае химического соединения характеристи¬ 
ческую температуру элемента Ѳ в атомарного вещества надо 
пересчитать на значение Ѳ с ' этого элемента в соединении по 
формуле Корефа 



где Т^ л и Т ва — абсолютные температуры плавления соеди¬ 
нения и элемента (хотя бы в обычных условиях и газообраз¬ 
ного). Зная характеристическую температуру Ѳ в ' и темпе¬ 
ратуру Т, при которой вычисляется теплоемкость, находят 
аргумент дебаевской функции Ѳп'/Т и но таблицам Дебая 
определяют теплоемкость с ѵ для каждого из сортов атомов, 
входящих в соединение (табл. 7). ' 

Характеристические температуры многих элементов при¬ 
водятся в справочниках. Если такие данные отсутствуют, ха¬ 
рактеристическую температуру элемента вычисляют по фор- і 
муле Линдемана: 



где Т пл — температура плавления; А — атомный вес; Е пл — 
атомный объем при температуре плавления. 

Атомный объем при температуре плавления рассчиты¬ 
вают, исходя из значения этой величины при комнатной 
температуре 

1 / пл= і/ 298 + МГ-298), 

где а ѵ — коэффициент теплового расширения. 

При отсутствии данных по коэффициенту теплового рас¬ 
ширения в формулу Линдемана подставляют атомный объем 
при комнатной температуре. Из-за малой величины коэффи¬ 
циента теплового расширения это, как правило, не вносит 
существенной ошибки в расчет характеристической темпера- 
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Таблица 7 


Значения дебаевской теплоемкости 


Ч Т 

щш 

Чіт 

шщж 

ч т 

с у* 

кал 

моль«град 

моль-град 

моль-град 


5,96 

3,2 

3,74 

7,8 



5,96 

3,3 

3,65 

8,0 

0.831 

0,2 

5,95 

3,4 

3,55 

8,2 

0,791 

0,3 

5,94 

3,5 

3,45 

8,4 

0,751 

0,4 

5,92 

3,6 

3,36 

8,6 

0,713 

0,5 

5,89 

3,7 

3,27 

8,8 

0,675 

0,6 

5,86 

3,8 

3 18 

9,0 

0,637 

0,7 

5,82 

3,9 

з;оэ 

9,2 

0,599 

0,8 

5,78 

4,0 

3,00 

9,4 

0,561 

0,9 

5,74 

4,1 

2,91 

9,6 

0,524 

1.0 

5,68 

4,2 

2,83 

9,8 

0,487 

1.1 

5,62 

4,3 

2,75 

10 

0,451 

1,2 

5,55 

4,4 

2,67 

11 

0,343 

1,3 

5,49 

4,5 

2,59 

12 

0,267 

1,4 

5,42 

4.6 

2,51 

13 


1,5 

5,35 

4,7 

2,43 

14 

0,168 

1.6 

5,27 

4,8 

2,35 

15 

0,137 

1,7 

5,19 

4.9 

2,27 

16 

0,113 

1,8 

5,10 

5,0 

2,20 

17 

0,094 

1,9 

5,01 

5,2 

2,06 

18 

0,079 

2,0 

4,92 

5,4 

1,94 

19 

0,067 

2,1 

4,83 

5,6 

1,82 

20 

0,058 

2,2 

4,74 

5,8 

1,70 

21 


2,3 

4,65 

6,0 

1,58 

22 

0,0433 

2,4 

4,55 

6,2 

1,48 

23 


2,5 

4,44 

6,4 

1,38 

24 

0,0335 

2,6 

4,34 

6,6 

1,30 

25 

0,0296 

2,7 

4,24 

6,8 

1,22 

26 

0,0262 

2,8 

.4,14 

7,0 

1,138 

27 

0,0234 

2,9 

4,04 

7,2 


28 

0,0211 

3,0 

3,945 

7,4 


29 

0,0189 

3.1 

3,84 

7,6 


30 

0,0172 


туры. Если справочные данные по атомному объему отсут¬ 
ствуют, его можно рассчитать из плотности твердого тела 


^298 — ’ 


Р298 


При рассмотрении спектра колебаний кристаллической 
решетки было показано (стр. 52), что характеристическая 
частота твердого тела ѵ т связана со скоростью распіюстра- 
нения колебаний в кристаллической решетке. Из этой фор¬ 
мулы следует, что 

'*т у 4гс 


■V. 
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Скорость распространения упругих колебаний в кристалли¬ 
ческой решетке совпадает со скоростью звука в твердом теле, 
поэтому характеристическая частота может быть рассчитана 
по скорости распространения звука ѵ в рассматриваемом 
веществе. Если такие данные отсутствуют, то скорость рас¬ 
пространения упругих колебаний можно вычислить из упру¬ 
гих констант кристалла. 

Известно, что в простейшем случае упругая волна в твер¬ 
дом теле — это совокупность трех волн: продольной и двух 
поперечных 

3 1,2 

Я I я 


В случае изотропных тел скорость распространения продоль¬ 


ной волны 






а поперечной волны 


В этих формулах р — плотность кристалла; (3 — коэффициент 
всестороннего сжатия; О — модуль сдвига. 

Теоретическому расчету поддается теплоемкость при по¬ 
стоянном объеме Су, для практических же расчетов необхо¬ 
димы данные по теплоемкости при постоянном давлении с р . 
Применение точной термодинамической формулы для пере¬ 
хода ОТ Су к с р 

Ср = с ѵ +ТѴ^- (20) 

требует знания коэффициента теплового расширения ау и 
коэффициента сжимаемости (3. Очень часто эти характери¬ 
стики твердых тел неизвестны, вследствие чего переход 
от Су к с р осуществляют, используя полуэмпирические зави¬ 
симости. Первая из этих зависимостей выводится из следую¬ 
щих соображений. По эмпирическому правилу Грюнайзена 
линейный коэффициент расширения твердого тела имеет ту 
же температурную зависимость, что и с р , т. е. 


где хі — некоторая постоянная величина. 
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Коэффициент теплового расширения (объемный) а ѵ свя¬ 
зан с линейным коэффициентом расширения соотношением 

«у=3а г . ( 22 ) 

Нернст и Линдеман нашли, что 



где Х 2 — постоянная величина. 

Подставляя выражения (21), (22) и (23) в формулу (20), 
получим уравнение Нернста — Линдемана для перехода от 
Су к с р через единственную характеристику вещества — тем¬ 
пературу плавления: 

С, = С„ + 0,0214-^- с/. (24) 

у 1 ПЛ 

Чаще всего уравнение Нернста — Линдемана решают ме¬ 
тодом последовательных приближений. Сначала подставляют 
в правую часть вместо с р рассчитанное по теории Дебая зна¬ 
чение Су и вычисляют с р ' в первом приближении. Далее под¬ 
ставляют вместо с р ее значение в первом приближении и вы¬ 
числяют уточненное значение с р " (второе приближение). 
Расчет повторяют до тех пор, пока значение с р , полученное 
из уравнения, не совпадает со значением из предыдущего 
приближения. 

Формулу Нернста — Линдемана можно представить в виде 
квадратного уравнения типа 

' ах 2 + Ьх + с = 0, 

где а = 0,0214 Т/Т п л ; Ь = — 1; с = с ѵ ; х = с р . 

Из теории квадратных уравнений известно, что уравнение 
такого типа не имеет вещественных корней, если 

Ь 2 — 4ас < 0, т. е. Ь 2 < 4ас. 

Для уравнения Нернста — Линдемана это условие прини¬ 
мает вид 

4-0,0214 ^ Гпл ѵ 

или 

-Т~с ѵ > 11 , 68 . 

В этом случае не имеет смысла решать уравнение и методом 
последовательных приближений, так как ряд расходится. 

Анализ полученного условия показывает, что пересчет 
изохорной теплоемкости на изобарную во всей области суще- 
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ствования твердой фазы, т. е. до температуры плавления, 
возможен только для одно- и двухатомных кристаллических 
веществ. В случае твердых многоатомных соединений об¬ 
ласть применимости формулы Нернста — Линдемана огра¬ 
ничена сверху предельной температурой, зависящей как от 
температуры плавления кристалла, так и от величины изо¬ 
хорной теплоемкости. Так как молярная теплоемкость соеди¬ 
нения равна сумме атомных теплоемкостей, а предельное 
значение атомной теплоемкости твердого тела равно 
6 кал/г* атом-град, то легко получить следующие темпера¬ 
турные границы применимости формулы Нернста — Линде¬ 
мана: 


Вещество 

Граница < 

применимости 

Одноатомное 

^ПЛ 

Двухатомное 

0,963 Т пл 

Трехатомное 

0,643 Т п л 

Четырехатомное 

0,483 Т пл 


Если с ѵ получено для температуры, превышающей ука¬ 
занный для данного вещества предел, ее пересчет на с р по 
формуле Нернста — Линдемана невозможен. 

Магнус и Линдеман нашли, что 



где хз — постоянная величина. Тогда 

С Р = С Ѵ + ТѴнѴ'Т. 

Объединяя постоянные величины в одну, получаем формулу 
Магнуса — Линдемана (или закон Т в степени три вторых) 

с р = с ѵ + аТ г ' 2 . 

Для того чтобы воспользоваться последней формулой для 
пересчета изохорной теплоемкости на изобарную, необходимо 
знать величину с р данного вещества при какой-нибудь одной 
температуре. Обычно известно значение изобарной теплоем¬ 
кости при 298 К (для многих веществ такие данные приво¬ 
дятся в і[7]). Поскольку Су при этой температуре легко рас¬ 
считать по теории Дебая, величину коэффициента а опреде¬ 
ляют из последней формулы. 

При полном отсутствии данных об изобарной теплоемко¬ 
сти величину коэффициента а можно оценить следующим 
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образом. Поскольку при небольших температурах (например 
при 298 К) разность между су и с р невелика, формулу Нерн¬ 
ста— Линдемана в первом приближении можно записать так: 

с р = с ѵ -\- 0,0214. 

и ѵ 1 ПЛ 

С другой стороны, по формуле Магнуса — Линдемана 

с р = с ѵ + аГЗ ' 2 - 

Приравнивая правые части двух последних уравнений, полу¬ 
чаем 


После того как определена величина коэффициента а, 
пересчет су на с р по формуле Магнуса — Линдемана может 
быть произведен для любой температуры. Значения функции 
Р /2 даются в табл. 8. 


Таблица 8 

Значения функции Т 3/2 


(к пересчету изохорной теплоемкости на изобарную 
по формуле Магнуса — Линдемана) 


Т 

7-3/2 

Г 

7-3/2 

298 

5 144,27 

1700 

70093 

300 

5196,15 

1800 

76 367 

400 

8 000,00 

1900 

82 819 

500 

И 180 

2000 

89443 

600 

14 697 

2100 

96234 

700 

18 520 

2200 

103 189 

800 

22 627 

2300 

110304 

900 

27 000 

2400 

117 575 

1000 

31 623 

2500 

125000 

1100 

36483 

2600 

132 574 

1200 

41 569 

2700 

140 296 

1300 

46872 

2800 

148162 

1400 

52 383 

2900 

156170 

1500 

58 095 

3000 

164 317 

1600 

64 000 




Для металлов разность между с р и су ничтожна при низ¬ 
ких температурах. При комнатной температуре эта разность 
составляет всего несколько процентов от величины теплоем¬ 
кости. Однако ее непременно надо учитывать при более высо¬ 
ких температурах. 
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Полупроводниковые материалы (элементарные полупро¬ 
водники и полупроводниковые соединения) характеризуются 
очень малым изменением объема с температурой, т. е. 
ничтожным значением коэффициента теплового расширения. 
Поэтому, как видно из табл. 9, для полупроводниковых ве¬ 
ществ с р практически не отличается от с ѵ . 


Таблица 9 


Теплоемкости полупроводниковых веществ 


Теплоемкость, 

кал/моль-град 

Алмаз 

81 

ОаАз 

ОаЗЬ 

Іп5Ь 

с ѵ 

1,44 

4,792 

5,50 

5,78 

5,90 

Ср 

1,46 

4,796 

5,52 

5,80 

5,91 

( с р —с ѵ ) эксперим. 

0,02 

0,004 

0,02 

0,02 

0,01 

(с р — с к ) по форму¬ 
ле (24) 


0,087 

0,13 

0,22 

0,28 


Таким образом, при расчете теплоемкости полупроводни¬ 
ковых материалов нельзя пользоваться формулой Нернста — 
Линдемана. Для них можно либо вообще пренебречь разни¬ 
цей между с р и с ѵ , либо пересчет с ѵ на с р осуществлять по 
формуле Магнуса — Линдемана. 

Последние теоретические исследования в области тепло¬ 
емкости привели к выводу формулы для пересчета с ѵ на с р , 
которая обеспечивает большую точность, чем уравнение 
Нернста — Линдемана, и позволяет рассчитывать с р полу¬ 
проводниковых соединений, для которых формула Нернста — 
Линдемана непригодна. Точное термодинамическое урав¬ 
нение, связывающее с ѵ с с р , можно представить в виде 


с р /с ѵ =1 + ТѴ 


(За,)» 

і*<ѵ ' 


С помощью уравнений динамики упругих сред показано, что 

ѴЗо ч _ 1 

$с ѵ ~ а • 


Следовательно, 
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где а — плотность упаковки атомов в кристаллической ре¬ 
шетке, т. е. отношение объема решетки, не занятого атомами, 
к общему объему. 


Тип решетки 

Плотность упа¬ 
ковки а 

Объемноцентрированная ку- 


бическая . 

0,68 

Гранецентрированная кубиче- 


ческая . . 

0,74 

Гексагональная . 

0,74 

Алмаза. 

0,34 


Пример. Рассчитать теплоемкость интерметаллического 
соединения Сг5і 2 при 1000 К по теории Дебая. Температуры 
плавления хрома, кремния и дисилицида хрома соответствен¬ 
но равны 2173, 1700 и 1773 К. Плотности хрома и кремния 
при комнатной температуре соответственно равны 2,33 и 
7,0 г/см 3 , а молекулярные массы — 52,00 и 28,08 г. 

По формуле Нейманна — Коппа: 

<Ѵ(СгЗі 2 ) = МСг) + 2<ѵ(3і). 


Для расчета атомных теплоемкостей по теории Дебая не¬ 
обходимо знать характеристические температуры элементов. 
Рассчитаем их по формуле Линдемана, предварительно вы¬ 
числив атомные объемы хрома и кремния: 

Ѵсг — ~уо~ — 7,43 см 3 /г-атом; 

Кзі = ТЗЗ~ =12 ,04 см / г ’ ат0м - 



Пересчитаем эти значения на характеристические темпера¬ 
туры элементов в соединении Сг5і 2 по формуле Корефа: 


ѳ сг = 437/-Щ- = 403 К; 



466 К. 


5 
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Далее вычисляем аргументы функции Дебая и по таблице 
находим соответствующие значения дебаевской теплоемкости: 

Ѳ Сг 403 

-у - = щд =0,403, с ѵ = Ъ,\\5 кал/г-атом-град; 
в!; 466 

~т~ ~ 1000 = 0,466, с ѵ = 4,984 кал/г-атом-град. 


Таким образом, изохорная теплоемкость соединения Сг8Ь 
с ѵ = 5,92 2-5,90= 17,72_кал/моль • град. 

Пересчитаем ее на изобарную теплоемкость по формула 
Нернста — Линдемана 

с р = 17,72+ 0,0214—с/ 


или . • 

0,01207 с р 2 — с р '-\- 17,72 = 0. 

Решая полученное квадратное уравнение, находим, что 
с р — 25,679 (корень уравнения 57,171 отбрасываем как не 
отвечающий физическому смыслу нашей задачи). 

Можно произвести пересчет с ѵ на с р методом последова¬ 
тельных приближений. Тогда получим (табл. 10) с ѵ — 25,677 
(пересчет произведен на ЭВЦМ «Наири»), Однако трчный 
пересчет без применения машинной техники отнимает много 
времени. 

Таблица 10 


Пересчет с ѵ нас р для СгЗі 2 по формуле Нернста — Линдемана 
методом последовательных приближений 


Номер прибли¬ 
жения 

с р , кал/моль-град 

Номер прибли¬ 
жения 

Ср, кал/моль-град 

і 

21,5099 

30 

25,6330 

2 

23,3044 

11 

25,6505 

3 

24,2751 

12 

25,6614 

4 

24,8326 

13 

25,6681 

5 

25,1630 

14 

25,6723 

6 

25,3624 

15 

25,6749 

7 

25,4840 

16 

25,6765 

8 

25,5586 

17 

25,6775 

9 

25,6046 

18 

25.6775 


Таким образом, расчетное значение с р = 25,68 кал/моль X 
X град. Расхождение с экспериментальным значением 
(24,41 кал/моль • град) составляет 5,2%. 
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Замеченные опечатки 


Строка Следует читать 

1-я снизу д// о 298 ( 8ІС) газ^ = И2 + 4 . 29 _ 4 . 97 = 

= — 160 ккал/моль. 

5- я снизу “ — 3,26 кал/град. 

6- я сверху \с р = — 5,03 + 0,83-10~ 3 Т + 

+ 1,71-КРГ* 2 кал/град. 

8-я снизу олова Заі*, если известны теплоемкости 
следующих жидких... 





















